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"  igh  degree  of  accuracy  in  order  to  be  able  to  cancel  the  interfering 
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SUMMARY 


This  report  documents  the  results  of  a  study  of  techniques  for 
cancelling  nonlinearity  generated  interference  from  the  input  of  sensitive 
receivers.  Harmonic  interference  is  caused  when  a  fundamental  signal  from 
a  communications  transmitter  impinges  on  a  nonlinearity  of  some  type, 
producing  harmonics  which  couple  into  the  antenna  of  a  colocated  receiver 
tuned  at  a  different  frequency.  More  complex  forms  of  interference  are 
produced  when  two  or  more  fundamental  signals  mix  in  a  nonlinearity, 
generating  intermodulation  products  which  couple  into  the  receiver.  .Means 
of  cancelling  both  forms  of  interference  by  utilizing  only  the  fundamentals 
have  been  investigated. 

The  study  consisted  of  an  analysis  of  the  interference  generated 
by  a  nonlinearity,  laboratory  experiment  to  verify  the  assumption  of  a 
memoryless  nonlinearity  and  to  evaluate  the  basic  approach  to  cancellation, 
and  a  computer  simulation  of  an  harmonic  interference  cancellation 
technique.  Various  cancellation  system  configurations  were  developed  for 
the  cancellation  of  harmonic  and  intermodulation  interferences  produced  by 
multiple  ncnlinearities.  The  laboratory  experiments  and  computer 
simulations  indicate  that  the  harmonic  cancellation  technique  appears  to  be 
realizable  and  practical.  The  report  concludes  with  a  discussion  of  some 
of  the  major  issues  that  must  be  resolved  in  any  practical  realization. 


The  author  acknowledges  the  assistance  of  burton  S.  Abrams  for 
some  of  the  analytical  and  experimental  results  contained  in  the  report. 
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EVALUATION 


Harmonics  and  inter  modulation  products  (IMP)  from  high  power  HF/VHF/UHF 
transmitters  aboard  B-52,  FB- 111,  and  CJ1  aircraft  are  interfering  with  the 
reception  of  a  satellite  starved  desired  UHF  signal.  (For  example,  the  tenth 
harmonics  of  an  HF  transmitter  (25  MHz)  falls  in  the  UHi-  range  (250  MHz).  Since 
the  interfering  and  desired  signal  occupy  the  same  frequency  band,  conventional 
filters  will  not  solve  the  problem  nor  will  conventional  Interference  Cancellation 
Systems  (ICS's)  solve  this  problem.  Conventional  collocated  ICS's  operate  on  the 
principle  that  a  small  portion  of  the  collocated  interfering  transmitter  signal  is 
sampled  and  adjusted  in  amplitude  and  phase  so  that  it  cancels  out  the  interfering 
signal  arriving  at  the  receiver  input.  However,  in  this  case,  the  "harmonic"  or 
"IMP"  interference  waveform  arriving  at  the  receiver  is  nonlinearly  related  to  the 
one  or  the  other,  or  both  offending  transmitters.  For  this  case,  a  harmonic 
cancellation  system  is  required  for  which  no  other  solution  is  apparent. 

The  significance  of  this  effort  is  that  it  has  direct  application  to  fighters, 
bombers,  and  C  I  aircraft,  reference  TPR  RAC.  In  addition,  this  technology  will 
have  direct  application  to  any  future  system  wherein  jammers  and/or  satellites 
common icat  ons  are  employed. 

h.  cM  a  n  z  V 

WAYNE  E.  WOODWARD 
Project  Engineer 
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1.0  i.\TR0i>UCT10.\ 

The  continually  inu'caair.j,  communication  requirements  aboard  Air 
force  airborne  platforms  of  all  types  is  crea l infa  a  continually  more 
critical  electromagnetic  compatibility  problem.  The  replacement  of  older 
communications  equipment  with  new  equipment  utilizing  solid  state 
components,  electronically  tuned  filters,  and  wider  band  coverage  is 
further  complicating  the  compatibility  problem.  Higher  transmit  powers  and 
new  transmission  waveforms,  especially  spread  spectrum  types,  add  to  the 
problem.  The  increasing  use  of  satellite-cased  communication  links  is 
surfacing  new  interference  problems  because  these  linscs  typically  operate 
with  little  margin. 

The  interference  problems  can  generally  be  divided  into  two 
classes:  linear  and  nonlinear.  Linear  interference  problems  are  those 
where  the  fundamental  signal  from  a  transmitter  couples  into  a  receiver  and 
prevents  the  receiver  from  detecting  or  demodulating  its  link  signal  with 
the  required  fidelity  or  sensitivity.  Nonlinear  interference  problems  are 
those  where  the  interference  coupled  into  the  receiver  is  a  harmonic  of  one 
transmitter  or  some  intermodulation  mix  of  two  or  more  transmitted  signals. 
The  nonlinear  interference  is  produced  by  one  or  more  transmit  signals 
impinging  on  a  nonlinearity  which  is  assumed  to  be  outside  the  victim 
receiver  itself. 

The  linear  interference  problem  can  De  solved  by  a  variety  of 
techniques.  Standard  tiitering  approaches  are  by  far  the  most  widely  used 
means  of  eliminating  this  form  of  interference.  in  addition,  active 
interference  cancellation  techniques  can  be  applied  to  solve  those  linear 
interference  problems  which  cannot  be  solved  by  the  application  of  spectral 
filtering.  While  not  as  mature  as  spectral  filtering,  linear  interference 
techniques  are  well-known  and  the  technology  is  available  to  solve  many 
problems.  Current  activity  in  linear  interference  cancelling  system 
technology  is  concerned  with  cancelling  higher  power  interference,  with 
broader  bandwidths,  taster  response,  less  complexity,  etc. 

Currently  available  techniques  of  solving  nonlinear  interference 
problems  are  limited  to  spectrum  management  and  the  elimination  of  the 
interference  producing  nonlinearities.  both  of  these  approaches  are 
becoming  increasingly  difficult  to  apply  because  of  the  demand  for  more 
communication  channels  and  the  wider  bandwidth  of  the  channels  as  well  as 
the  trend  towards  broader  band,  electronically  tuned  communications 
equipment  which  frequently  is  the  source  of  the  nonlinearity.  In  order  to 
provide  an  alternate  technique  for  combatting  nonlinear  interference,  the 
application  of  active  cancellation  concepts  to  nonlinear  interference 
problems  Is  being  studied. 

This  report  documents  the  results  of  a  study  of  the  nature  ot 
nonlinear  interference  from  the  standpoint  of  applying  active  cancellation 
techniques  to  its  reduction.  The  study  involved  analytical  work, 
experimental  measurements,  and  computer  simulations  to  characterize  the 
problem  and  to  evaluate  techniques  for  the  cancellation  of 
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nonlinearity-generated  interference.  Most  of  the  study  was  directed  at 
nonlinear  interference  which  consists  of  harmonics  of  a  single  fundamental, 
since  this  is  less  complex  than  the  more  general  problems  involving  mixes 
of  several  transmitter  waveforms.  An  understanding  and  approach  to  solving 
harmonic-type  interference  will  lead  to  a  solution  of  the  more  general 
nonlinear  interference  problems. 

An  approach  to  providing  cancellation  of  harmonic  interference 
lias  been  developed  and  evaluated  by  computer  simulation.  This  approach 
involves  the  synthesis  of  the  envelope  of  the  harmonic  interference  by 
processing  the  fundamental  envelope  and  the  received  harmonic  envelope. 
The  synthesized  envelope  is  used  to  modulate  a  carrier  at  the  harmonic 
frequency  which  is  then  used  to  effect  cancellation  of  the  received 
harmonic  interference. 

The  canceller  system  is  implemented  by  utilizing  standard 
canceller  techniques  plus  harmonic  synthesis  processing  consisting  of 
straightforward  digital  and  analog  signal  processing  component-.  It  is 
recommended  that  an  experimental  model  of  this  equipment  be  developed  in 
order  to  further  evaluate  the  capabilities  of  such  techniques. 

1.1  BACKGROUND 


The  identification,  evaluation  and  elimination  of  nonlinearities 
on  board  Air  Force  aircraft  is  addressed  in  detail  in  reference  [1].  As 
described  in  i 1 J ,  a  series  of  measurements  were  made  on  several  aircraft. 
The  measurements  indicated  that  the  predominant  harmonic  interferences  were 
not  produced  by  airframe  effects  ("rusty  bolts,"  mechanical  joints,  etc.) 
as  initially  expected,  but  by  nonlinearities  in  the  communication  and 
navigation  systems. 

Further  tests  and  measurements  located  some  of  the  most 
troublesome  causes  of  these  nonlinearities  and  the  relative  magnitude  of 
the  distortion  produced  by  these  sources.  Equipment  modifications  to 
reduce  the  production  of  distortion  products  were  investigated, 
implemented,  and  evaluated. 

Subsequent  measurements  of  the  AFSATCOM  receiver  have  indicated 
that  additional  nonlinearities  are  present  which  produce  distortion 
products  which  fall  in  the  satellite  downlink  band,  thereby  limiting  the 
sensitivity  of  the  receiver.  This  loss  of  sensitivity  is  particularly 
critical  for  satellite  links  such  as  AFSATCOM  or  the  Global  Positioning 
System  (NAVSTAR)  since  these  links  are  usually  operated  with  little  margin. 
Aircraft  on  which  distortion  product  interference  is  causing  problems  are 
the  EC-135C,  RC-135,  FB-lllA,  and  B-52U.  References  [2],  [3],  and  [4] 
contain  descriptions  of  other  studies  which  analyzed  nonlinearities  and 
means  of  locating  them  on  structures  as  large  and  complicated  as  Navy 
ships.  The  general  problem  of  analysis  of  nonlinear  effects  has  been 

[1]  "Nonlinear  Interference  Cancellation  System,"  RADC-TR-78-225 ,  November 
1978. 
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address  d  by  researchers  in  various  fields  for  many  years.  The  early  work 
of  Wiener  as  well  as  tuc  ..it  -r  work  t5J  represents  soi.ie  «£  the  earliest 
efforts  to  apply  th-->  tools  of  mathematical  analysis  to  the  solution  uL 
nonlinearity  problems  in  engineering.  Much  of  the  analysis  can  be  divined 
into  consideration  of  two  basic  types  or  nonlinear  systems:  those  with  no 
memory  and  those  containing  memory.  Memory less  systems  (both  linear  .mb 
nonlinear)  have  the  property  that  the  system  output  .it  each  instant  ot  time 
is  determined  by  the  input  at  that  instant  and  is  unaffected  by  the  past 
inputs  to  the  system. 

Much  of  the  earlier  applied  work  was  restricted  to  tiu  analysis 
of  memoryless  types  of  nonlinearities.  More  recent  efforts  nave  extended 
the  earlier  work  to  include  nonlinearities  with  nePiory.  Analysis  of  these 
more  complicated  nonlinearities  requires  more  powerful  analytical  tools. 
Recently,  there  has  been  renewed  interest  in  the  solution  of  these  problems 
by  application  of  Volterra  functional  series.  Reference  h  contains  a 
description  of  this  series  and  its  application  Lo  nonlinear  analysis  and  is 
typical  of  current  applications  of  Volterra  analysis. 

Most  work  reported  to  date  involves  the  analysis  of  the 
distortion  characteristics  of  devices  ana  the  computat ; -n  of  the  levels  of 
distortion  products  (UP)  produced  by  tiiese  devices  or  the  exact  fundamental 
wave  for  ss  in  the  nonlinearity.  Lsually,  tiie  analysis  ot  the  distortion 
chur.ict  ristics  of  these  devices  is  conducted  for  CW  signals  and, 
therefore,  is  not  direct Ly  applicable  to  the  problem  of  a  more  complex 
nodulated  carrier.  No  work  is  reported  that  applies  directly  to  the 
problem  of  cancelling  UP's  at  harmonic  or  intermodule tion  products  jjroducee 
by  nonlinear  devices  excited  by  one  or  two  incident  signals. 


L  2  J  Riga,  Walter  H.,  "Spurious  Signals  Generated  by  electron  Tunneling  on 
Large  Reflector  Antennas,"  Proc  IEEE,  vol.  !  \  no.  2,  February  197b. 

[3]  Chase,  W.M. ,  J.W.  Rockway  and  G.C.  Salisbury,  "A  Method  of  Detecting 
Significant  Sources  of  Intermodulation  Interference,"  IEEE  Trans  or. 
Electromagnetic  Compatibility,  vol.  EMC-17,  no.  2,  May  1975. 

[4J  Chase,  W.M. ,  "Ship  RF I  Survey  Procedure  for  liF  Frequencies,"  NELC 
Technical  Document  336,  21  June  1974. 

[5]  Wiener,  N.  Nonlinear  Problems  in  Random  Theory,  MIT  Press,  1958. 

(6j  "Nonlinearity  System  Modelling  and  Analysis  with  Application  to 
Communication  Receivers,"  kADC  TR-73-178,  June  1978. 

[7J  Thomas,  E.J. ,  An  Adaptive  Echo  Canceller  in  a  Nonideal  Environment 
(Nonlinear  or  Time  Variant),"  bS'l  J ,  vol.  5u,  no.  8,  October  1971. 

[8]’  Ijlconer,  D.D.,  Adaptive  Equalization  ot  Channel  Nun  1 i nea ri t ies  in 
n.Vi  Lata  Transmission  Systems,"  bSTJ ,  Vol.  57,  no.  7,  September 
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Some  work  is  underway  to  provide  equalization  for  the  nonlinear 
characteristics  of  telephone  channels.  References  [7]  and  [8]  report  on 
the  application  of  these  techniques  to  baseband  equalization  of  telephone 
line  nonlinearities.  The  work  reported  in  reference  [8]  compensated  for 
quadratic  and  cubic  nonlinearities 

(amplitude  only)  with  a  resulting  improvement  in  the  bit  error 
rate  performance  of  the  channel.  Reference  [7]  is  a  mathematical  analysis 
and  simulation  of  the  performance  of  an  equalizer  when  confronted  with 
nonlinear  distortion  produced  by  phase  jitter  caused  by  local  oscillator 
frequency  instability  in  the  transmission  channel. 

The  approaches  to  DP  cancellation  considered  during  this  study 
assume  that  the  nonlinearic ies  to  be  encountered  in  practice  are  memory less 
and  continuous  with  no  hysteresis  or  jump  characteristics.  Laboratory 
measurements  appear  to  support  the  truth  of  these  assumptions.  Typical 
communication  system  components  such  as  diodes,  amplifiers,  mixers  and 
electronically  tuned  radio  front  ends  have  been  used  in  the  experiment. 
The  measurements  made  indicate  that,  at  least  over  the  bandwidth  of 
interest  (tens  of  kilohertz),  the  nonlinearly  generated  harmonic  products 
are  independent  of  frequency,  implying  that  the  nonlinearities  have 
sufficiently  small  memory  effects  that  they  can  be  regarded  as  memoryless. 
It  is  Likely  that  as  larger  frequency  bands  are  considered  (large 
percentage),  memory  effects  would  become  more  significant. 

It  should  be  enpnasized  that  the  lack  of  memory  effects  described 
above  refers  to  narrow  bands  centered  around  a  particular  harmonic  or 
intermodulation  product  carrier  frequency  and  does  not  refer  to  wide  bands 
spanning  the  radio  frequency  difference  between  the  various  harmonics. 
Therefore,  mathematically,  these  nonlinearities  are  not  strictly 
memoryless.  This  distinction  will  be  addressed  again  later. 
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In  this  section  the  harmonic  and  in te modulation  product  signal;, 
generated  when  one  or  more  modulated  fundamental  carrier  signal**  are 
incident  on  a  nonlinearity  are  analytically  investigated.  In  addition,  the 
effects  of  multiple  nonlinearities  are  investigated.  As  mentioned  earlier, 
the  analysis  assumes  that  the  nonlinearities  are  memory  less ,  continuous  and 
single-valued.  This  assumption  appears  to  be  justified  based  on 
measurements  oi  the  nonlinearity  characterisL  ics  01  typical  communications 
equipment  and  components.  Section  3  contai ns  a  summars  oi  those 
measurements. 

The  distortion  product  analysis  is  divided  into  three  sections. 
First,  harmonics  of  a  single  fundamental  signal  impinging  on  a  nonlinearity 
Mt  considered.  Next,  the  iniormodulation  distortion  products  prouuced 
when  two  signals  combine  on  a  nonlinearity  are  analyzed.  Finally,  the 
effects  of  multiple  nonlinearities  on  the  characteristics  of  Lite  total  bP 
interference  encountered  by  a  receiver  are  analyzed. 

2.1  HARMONIC  DISTORTION 

Meddling  non  1  incur  it  ics  as  continuous,  single-valued  and 
memory’..  *  .  i  lows  the  MLj-ui  of  u  itonl  iiiearity  to  be  expressed  .is  a  power 

series  ex,v.i  iron  or  the  incident  fundamental  signal.  The  output  ot  such  a 
nonlinearity  is  given  by 

OO 

v0(t)  =  l  ak  Vl(t)  (2-1) 

k=0 


where  v  (t)  is  the  incident  fundamental  signal.  The  characteristics  of  a 
particular  device  are  contained  in  the  set  of  coefficients  a^.  Tne 
memoryless  nature  of  a  device  characterized  by  an  expression  such  as  (2-1) 
is  apparent  since  the  output  at  any  time  depends  only  on  the  input  at  that 
time  and  is  independent  of  the  past  history  of  the  input.  For  the  same 
reason,  a  device  whose  behavior  is  described  by  (2-1)  cannot  display  any 
hysteresis. 


The  fundamental  signal  incident  on  the  nonlinearity  will  be 
represented  by 

v,(t)  **  V(t)  cos(u)Qt+0(t))  (2-2) 

This  representation  includes  amplitude  and  phase  modulated  carriers  of 
various  types,  depending  on  the  selection  of  the  amplitude  function  V(t) 
and  the  ph. >.*,«.  tunction  j(t).  Inserting  (2-2)  in  (2-1)  yields 

vO 

vnU)  =  L  ai,  v’k(t)  cosk(w  ( t ) )  (2-3) 

k=0  u 
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Terms  in  this  sum  may  be  expanded  using  the  results  in  {. 9 J  to  give 


where 


with 


V  (t)cos  ((jo^t+iji)  =  v  (t)  £  2C(k,m)cos  (itta.t+m<J>) 

u  m=0  u 

m  even  if  k  even 
m  odd  if  k  odd 


C(k,m)  = 


e  r(k+D 

tn 


2k+1r(i  -  ~>r(i  + 


£  =1  for  m  =  0 

m 

=  2  for  m>0 
f(k)  -  Gamma  function 


=  (k-1) !  for  integer,  k,  k>0. 


(2-4) 


by  utilizing  che  properties  of  the  Gamma  function,  it  can  be  shown  thaL 
c(k,m)  =  0  for  n>k. 

Using  (2-4;  in  (2-3)  allows  the  output  of  the  nonlinearity  to  be  expressed 
as 


v0(t)  =  2  J  a.V  (t)  £  C(k,m)cos (moj^t+m^)  (2-5) 

k=0  m=0 

(k-m)  even 

or  by  interchanging  summations  as 

00  DO 

vQ(t)  =  2  l  cos(ntk>0t+m<f>)  l  akVk(t)C(k,m)  (2-6) 

m=0  k=m 

(k-m)  even 

Uadi  of  the  above  representations  displays  the  structure  of  the 
harmonic  distortion  products  in  a  different  way  Equation  (2-5)  shows  that 
each  term  in  the  power  series  expansion  of  the  nonlinearity  develops  a  sum 
of  carriers  at  (alternate)  harmonics  of  the  fundamental  frequency  and  that 
the  fundamental  phase  term  <Mt)  is  scaled  by  the  harmonic  number.  All 
these  terms  are  amplitude  modulated  by  VK(t).  Kor  example,  for  the  k-9 
tern,  the  following  components  are  present: 


1 9 j  Davenport ,  W. B. ,  Root,  W. L. ,  An  Introduction  to  t he  Theory  of  Random 
■Signals  and  Noise ,  McGraw-Hill  book  Co.,  New  York,  1978,  p.  284. 


V9(t){c(9,l)cos(;j0t+<‘>)  +  C(9,3)cos(3a)yt+3<{))  +  C(9,5)cos(5(jj^t+54>) 
+C(9,7)cos(7u)Qt+7(J))  +  C(9,9)r.os(9ijQt+9^)  }  (2-7) 


This  representation  shows  which  harmonic  carriers  are  amplitude 
modulated  by  \r(t).  The  alternate  representation  of  (2-6),  on  the  other 
hand,  identifies  the  signal  structure  at  each  harmonic.  This  can  be  seer, 
by  considering  the  second  sum  in  (2-6)  for  a  particular  m.  This  sum,  which 
is  a  combination  of  powers  of  V(t),  amplitude  modulates  the  n-th  harmonic 
carrier.  The  bandwidth  of  the  modulation  of  the  r.i-th  harmonic  it, 
detrermined  by  the  bandwidth  of  the  fundamental  modulation  V(t)  and 
resulting  bandwidth  of  the  powers  of  V(t)  in  the  sum.  The  energy  in  those 
higher  power  terms  diminishes  rapidly  as  k  increases  because  of  the  big), 
order.  Thus  the  bandwidth  of  the  m-th  harmonic  energy  is  concentrated 
about  the  m-th  harmonic  carrier  frequency.  For  the  fundamental  carrier 
frequencies  under  consideration  here  (2  Mhz  or  higher),  there  typically  is 
no  overlap  of  energy  between  adjacent  harmonics. 

Equation  (2-b)  clearly  shows  the  structure  of  these  harmonics. 
Generally  a  narrowband  receiver  is  tuned  to  receive  only  one  of  these 
harmc  •cv.  The  harmonic  interference  at  the  n-th  harmonic  is  given  by 

vn(t)  =  Vn(t)  cos (no)gt+n<j> ( t ) )  (2-8) 

where  the  amplitude  modulation  vn(t)  is 

CO 

Vn(t>  =  2  l  akC(k,n)Vk(t) 
k=n 

(k-n)  even 

As  (2-8)  shows,  the  amplitude  modulation  of  the  n-th  harmonic  carrier  is  a 
sum  of  alternate  powers  of  the  fundamental  modulation  V(t),  beginning  with 
the  n-th  power.  The  terms  in  the  sum  are  weighted  by  the  nonlinearity 
coefficients  a^  and  the  expansion  coefficients  C(k,n).  As  described 
above,  the  bandwidth  of  the  sidebands  at  the  n-th  harmonic  is  determined  by 
the  fundamental  modulation  spectrum  and  the  coefficients  a.  and  c(k,n). 

In  principle,  this  bandwidth  is  infinite,  as  shown  by  theK  expansion  in 

(2-8)  tor  V  (t).  However,  the  energy  near  the  band  edges  of  Vk(t) 
diminishes  rapidly  for  large  k  since  the  spectrum  of  Vk(t)  is  the  k-fold 
convolution  of  the  spectrum  of  V(t).  In  general,  the  dominant  term  in  the 
sum  is  the  k=n  term,  since  higher  powers  of  V(t)  will  be  smaller  in 
ma;;n  i  r.  ode  . 
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2.2  INTERMODULATION  PRODUCT  DISTORTION 

In  this  section  the  analysis  of  the  distortion  products  produced 
by  a  memory less  nonlinearity  is  extended  to  the  case  where  two  amplitude 
modulated  fundamentals  are  incident  on  the  nonlinearity.  The  output  of  the 
nonlinearity  is  given  by  (2-1)  where  now  the  input  is  given  by 

v(t)  =  A^  (t)cosu^t  +  (tjcos^t  (2-9) 

in  which  A.(t)  and  A£(t)  are  amplitude  modulations  for  carriers  at 
frequencies  f^  and  fo.  For  simplicity  only  amplitude  modulation  is 
considered.  Inserting  T2-9)  in  (2-1)  yields 

00 

vQ(t)  =  I  an[A1coscJ1t+A2cosa)2t]n 
n=0 


oo  n 

=  1  Bn  1  ^i^AiA2  icosi(i)1tcosn~i(i)2t  (2-10) 

n=0  i=0 

where  the  dependence  of  A.  and  A,  on  t  has  been  dropped  for  convenience 
and 

(k)  =  _ 111 _ 

V  i ! (k-i) ! 


are  the  binomial  coefficients. 


In  order  to  show  the  intermodulation  frequencies  explicitly, 
(2-10)  will  be  expanded  further  by  use  of  the  following  [10]: 

i  n 

cosn0  =  —  l  (?)  cos(n-2j)0 
2n  j=0  ^ 

Using  this  in  (2-10)  results  in 


[10]  Jolley,  B.W.,  Summation  of  Series ,  Dover  Press,  1961. 
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v(J(t) 


I* 

i 

1 

C 

a 

T 
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r 

n  ,n 

i=0 

j=0 

k=0 

2n  1 

n 

l 

i 

I 

n-i 

I 

a 

n 

i=0 

j=0 

o 

n 

2*1+1 

+  (n-i-2k)(jj2t]+cos[  (  i-2  j  (r>-i-2k)w0t  ]  } 


(2-11) 


which  shows  that  the  in termoou ration  products  generated  by  the  nonlinearity 
occur  at  frequencies 

f  j  =  r.i  j  f  j  +  m9f  2 
where  from  (2-11) 

m!  =  i~2j 

i«2  =  n-i-2k 

In  general,  re  and  range  over  the  set  of  positive  anti  negative 
integer**  c*--  *un  be  seen  by  examination  of  the  limits  in  (2-11).  The  order 
of  the  i  i.* >  r.  *>dulat  ion  product  (IMP)  is 

IMP  order  =  |mj  +  |m2! 

Since  we  are  considering  IMP  interference*  to  a  narrowband 
receiver  tuned  to  a  frequency  f  it  is  not  necessary  to  retain  all  the 
terms  in  (2-11).  Only  those  IMP  which  occur  at  f  will  get  through  the 
receiver  passband.  Furthermore,  if  we  assume  that  the  modulation  sidebanus 
of  the  IMP  fall  off  rapidly  so  that  there  is  no  appreciable  spectral 
overlap  betweeen  the  IMP,  only  those  IMP  occurring  at  f  will  arfect  the 
receiver.  Thus,  the  IMP  interference  at  t  is  given  by  those  terms  in 
(2-11)  for  which 


MM 


rl2r  2  • 


There  are  multiple  solutions  for  Equation  (2-13),  that  is,  sets  of 

^mli’m2p  for  i  =  M2,...  which  yield  the  IMP  at  f^.  However,  the 

IMP  which  causes  the  greatest  interference  are  those  where  the  order  of  tne 

IMP  is  low,  typically  3,  3,  or  7.  Therefore,  for  a  particular  fR  ana 

fixed  fj  and  *2,  the  values  of  mj  and  r.»2  in  (2-13)  which  cause  the 
interference  can  be  considered  as  unique  and  will  be  such  as  to  provide  a 

low  order  IMP.  Hence,  we  may  consider  m  and  m2  as  fixed  and  apply  the 
definitions  of  (2-12)  to  reduce  (2-11)  to  that  IMP  which  occurs  at  f 

Making  the  subsitution  m.  =  i+2j  and  bm  =  n-i-2k  eliminates 
the  sums  over  k  and  j  in  (2-11)  and  yields  the  IMP  at  f  as 


VR(t) 


00  n-m„  a 

l  I  2  -^7 <?><* 


i  4n-i 


n*m^+m2  i^  W ^-^2  *c-[W*W  ] 


even  even 


(2-14) 


where  only  the  r.i^fj+n,  f„  term  has  been  retained  from  (2 
transform ing  the  summation  indices  using 


-11).  By 


^2  =  n-i-m2 

(2-14J  becomes 

£ . +m.  £  +m2 

00  °°  a£  +ra  +£?+tn  (£.+m1+£2+m2) !  A? 

V|!<t)  ‘  .  ln  A  «,«,'+*, V-'l  - «7_S -  '“»["1“lt+*2“2tl 

V°  l2-0  2  1  1  2  2  (J)!(jl4.1)!0l<j^)! 

Uneven  £2even 


1'ne  IMP  at  f  can  be  expressed  as 


vR(t)  =  y(t)  cosffm^j+m^Jtl 


(2-15) 


00  00  fc.+m.  S.,+tn2 

y(t)  -I  l  c.  „  A.1  i(t)A/  Z(t) 

V°  V°  12 

£^  even  £2  even 

where  the  constants  cj^,^  include  the  expansion  coefficients  and  the 
nonlinearity  coetficients  ana  are  given  by 

aJl1-Hn1+il2+in2  (^1-hn1+£2+m2) ! 

c£  ,£,  =  £1+m1+£7+m,+l  ~T~  ST  £~  17 

1  1  2  1  1  1  (2i)!(24m1)!(^)!(2^+ m2)! 

2.3  MULTIPLE  NONLINEAR I TIES 

The  underlying  assumption  for  the  previous  developments  is  that 
the  nonlinearity  generating  the  interference  is  strictly  memoryless. 
Nonlinear  devices  that  exhibit  AH— PM  conversion  are  not  memory less,  for  it 
takes  memory  to  vary  the  zero  crossings  of  the  output  waveform  from  those 
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of  the  input.  Sucli  an  efiect  can  bo  produced  by  a  nonlinear  device  as 
simple  as  a  diode  wit;,  a  capacitor  across  the  junction,  so  that  the  R-C 
time  constant  depends  on  the  input  amplitude.  More  importantly,  as  will  be 
shown,  interference  produced  by  a  composite  nonlinearity  that  is  formed  by 
a  group  of  memoryless  nonlinearities  whose  waveforms  are  combined  with 
different  phase  shifts  (e.g.,  due  to  different  propagat ion  delays)  aiso 
exhibits  AM- PM  conversion. 

As  a  practical  matter,  nonlinearities  are  likely  to  occur  in 
multiples  due  to  the  presence  of  multiple  communication  systems  on  board  a 
single  aircratt.  In  addition,  it  is  likely  that  a  single  communications 
terminal  (including  antenna,  tuner,  one  or  more  receivers  or  transmitters) 
will  contain  several  nonlinearities.  Such  a  composite  nonlinearity 
consists  of  individual  nonlinearities  radiating  harmonic  energy  that  is 
received  by  an  antenna  on  the  aircraft  witli  slightly  different  propagation 
delays.  This  model  is  called  weakly  dispersive  because  it  is  assumed  that 
the  envelope  delays  are  insignificant,  but  the  effect  on  carrier  phase  is 
not.  Thus ,  a  sum  of  harmonic  waveforms  is  received  which  constitutes  the 
total  harmonic  interference  that  may  be  expressed  as 

v  (t)  =  £  A^(t)  cos  (nui  t+ni})(t  )+6  .  j  (2-16) 

i 

where  A  t.t)  is  the  harmonic  amplitude  modulation  associated  with  each 
nonlinearity.  Each  tern  in  the  sum  of  (2-16)  corresponds  to  an  harmonic 
interference  as  derived  earlier  and  expressed  in  (2-8).  The  summation 
variable  i  is  used  to  index  each  of  the  contributing  nonlinearities.  The 
phase  shift  q  representing  the  effect  of  propagation  delay,  is 
subscripted  by  the  summation  index  i  to  show  that,  in  general,  its  value  is 
different  for  each  value  of  i.  Since  the  nonlinearities  may  vary  greatly 
in  their  particular  characterist ics  and  in  the  amplitudes  of  their  input 
fundamental  waveforms,  the  envelopes  A.(t)  will  also  differ  for  different 
valuesofi.  1 

Thus,  the  total  n-th  harmonic  can  be  expressed  as 

vn  (i )  =  lZ  A.(t)cos0  |cos[no)  t+n<Kt)]-[y  A.  (t)  sin0  ]sin[nu>  t+n<J>(t)]  (2-17) 
i  i  1  1  c 

in  which  the  inphase  and  quadrature  components  are  evident.  This  can  be 
made  more  apparent  by  writing 

vn(t)  =  )  cos  [nojct+n$(  t )  ]+Aq(t ) sin  [nu)c  t+n<{>(t )  ]  (2-18) 

where 

Aj(t)  =  £  Ai(t)  cosB. 


Aq(t)  =  ~l  A1(t)  sinO. 
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Thus  it  can  be  seen  that  mu  It  i ;» i  e  tior.l  inearities  produce  AM  to  PM 
conversion  which  must  be  accommodated  in  any  system  attempting  to  effect 
harmonic  cancellation.  The  AM-to-PM  conversion  is  evident  by  rewriting 
vn(t)  once  more  as 

vn(t)  =  A(t)  cos  [nu)^t+n<}> ( t  )+0  ( t )  ]  (2-19) 


where  the  envelope  A(t)  is 


A(t)  =  Aj (t)  +  Ap(t)  h 


anu  the  phase  9v.t)  is 


-1  V0 

0(t)  =  tan  1{  . 

AQ(t) 


1 


/ 


EXPERIMENTAL  RESULTS 


3  U 

j.l  FREQUENCY  DEPENDENCE  UP  NONLINEARITIES 

The  analysis  of  nonlinearity-generated  interference  of  either  the 
liarmonic  or  intermodulation  product  type  is  greatly  simplified  if  the 
nonlinearity  in  question  is  memoryless.  Likewise  the  structure  and 
complexity  of  a  harmonic  distortion  cancellation  system  is  greatly  reduced 
if  the  nonlinearity  is  memoryless.  The  nonlinearity  models  used  in  Section 
2  assumed  that  there  are  no  memory  effects  present  in  the  nonlinearity  and, 
as  a  consequence  of  this  assumption,  the  entire  liarmonic  spectrum  displayed 
no  memory  effects  at  all.  In  practice  this  is  unlikely  to  be  true. 
However,  all  that  is  really  needed  is  that  memory  effects  be  negligibly 
small  over  a  frequency  band  commensurate  with  the  bandwidth  of  the 
modulation  sidebands  of  the  harmonic  of  interest.  Lack  of  memory  over  a 
band  of  this  extent  is  more  likely  to  be  found  with  real  nonlinearities  and 
implies  that  the  amplitude  of  the  modulation  sidebands  are  not  a  function 
of  the  fundamental  modulation  frequency. 

In  order  to  experimentally  evaluate  these  effects,  two 
nonlinearities  were  tested  in  the  laboratory  to  determine  if  their  harmonic 
responses  to  an  AM  waveform  were  dependent  on  the  frequency  of  the 
fundamental  modulation.  A  block  diagram  of  the  test  setup  is  given  in 
Figure  i . 


A  sinusoidal Ly  modulated  All  double-sideband,  suppressed  carrier 
waveform  was  synthesized  by  summing  two  sinusoids  at  equal  levels  with  a 
slight  frequency  difference  A,  with  both  frequencies  set  near  25  MHz.  The 
sum  was  amplified  and  lowpass  filtered  for  harmonic  rejection.  The 
resulting  waveform  was  then  applied  to  the  test  nonlinearity,  and  the  sixth 
harmonic  at  150  MHz  was  monitored  on  a  spectrum  analyzer.  Two  test 
nonlinearities  were  used:  an  KD100  Schottky  diode  and  the  antenna  port  of 
an  ARC-164  radio  with  its  power  off.  The  attenuator  preceding  the 
nonlinearity  was  adjusted  to  obtain  approximately  equal  sixth  harmonic 
power  levels  for  both  nonlinearities. 

Both  the  ARC-164  and  the  diode  sixth  harmonic  response  were 
observed  with  approximately  I  kHz  and  10  kHz  tone  spacing  (at  the 
fundamental  frequency).  Figure  2  is  a  photograph  of  the  spectrum  analyzer 
display  with  the  ARC-164  used  as  the  test  nonlinearity.  The  upper 
photograph  is  for  a  1  kHz  spacing  between  the  two  fundamental  tones  (i.e., 
25.000  MHz  and  25.001  Mhz)  while  the  lower  photograph  is  for  a  10  kHz 
spacing  between  the  two  fundamental  tones  (25.000  MHz  and  25.010  Mhz). 
Figure  3  is  for  the  same  conditions  with  the  analyzer  vertical  sensitivity 
improved  to  2  dB/div.  Comparison  of  the  upper  and  lower  analyzer 
photographs  shows  that  there  is  no  discernible  difference  between  the 
relative  sideband  levels,  indicating  that  the  harmonic  sideband  structure 
is  not  dependent  on  the  frequency  of  the  fundamental  modulation.  The 
measurement  technique,  instrument  accuracy  and  sensitivity  permits 
measurement  of  differences  in  sideband  levels  as  small  as  0.5  or  0.25  dB, 
indicating  that  the  ARC-164  front  end  (unpowered)  is  memoryless  over  a  band 
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FIGURE  1.  TEST  SETUP  TO  INVESTIGATE  FREQUENCY  DEPENDENCE  OF  NONLINEARITIES 


1  kHz  FUNDAMENTAL 
TONE  SPACING 


VERTICAL:  10  dB/DIV 
HORIZONTAL:  2  kHz./DIV 


10  kHz  FUNDAMENTAL 
TONE  SPACING 


VERTICAL:  10  dB/DIV 
HORIZONTAL:  20  kHz/DIV 


FIGURE  2 


ARC-164  TEST  NONLINEARITY 


10  kHz  FUNDAMENTAL 
TONE  SPACING 


VERTICAL:  2  dB/DIV 
HORIZONTA:  20  kHz/DIV 


FIGURE  3 


ARC- 164  TEST  NONLINEARITY  -  DETAIL 


L  t.  1 1..  ul  K.1  l  \  }i  ix.  i  i...  w  i  d  l'  Lo  . 1  L  i  v.  .. L  '  .  )  to  0.23  d  to  • 

i  i.-,ur..'_.-  4  and  3  .  n  ■/.  the  sar.:e  sequence  of  analyzer  photographs 
with  the  diode  used  as  tile  test  nonlinearity.  The  sane  comparison  can  no- 
made  showing  that  diode  is  also  memoryless  over  the  band  measured. 

By  comparing  Figure  2  to  Figure  4  it  can  be  seen  that  the 
sideband  structure  of  the  AKC-164  nonlinearity  and  the  diode  nonlinearity 
are  different.  For  example,  comparison  of  the  lower  photographs  in  each 
figure  indicates  that  the  ARC-164  sidebands  fall  off  more  rapidly  than 
those  of  the  diode.  This  difference  in  relative  sideband  structure  means 
that  the  sixth  harmonic  waveform  from  the  diode  could  not  be  used  to 
achieve  significant  cancellation  of  interference  produced  by  an  ARC-164 
front  end. 

Appendix  A  analyzes  the  cancellation  that  can  be  achieved  by 
using  the  harmonic  distortion  produced  by  one  nonlinearity  to  cancel  the 
harmonic  distortion  produced  by  a  second  nonlinearity. 

3.2  TIME  DOMAIN  SYNTHESIS  OF  HARMONIC  ENVELOPE 

In  order  to  evaluate  the  concept  of  harmonic  cancellation  by 
synthesis  of  the  harmonic  envelope,  a  laboratory  experiment  was  conducted 
using  s'.  -tndard  laboratory  equipment  and  an  interference  cancellation 
system.  t..>icully ,  die  experiment  consists  of  generating  a  sinusoidally 
amplitude  modulated  fundamental  which  excited  a  nonlinearity  to  generate  a 
third  harmonic  interference.  Direct  synthesis  of  the  envelope  of  a  third 
hjrmonic  signal  to  effect  cancellation  of  the  interference  was  used.  The 
synthesis  ot  the  third  harmonic  envelope  was  accomplished  by  use  of  a 
Wavetek  waveform  generator,  the  output  of  which  was  used  to  amplitude 
modulate  a  carrier  at  (tie  third  harmonic.  Cancellation  was  acccompli shed 
by  use  of  the  MX-500,  a  UI1F  cancellation  system  with  adaptive  control  loops 
which  is  capable  of  providing  60  db  of  cancellation. 

The  waveform  synthesizer  used  was  a  Wavetek  Model  No.  175, 
Arbitrary  Waveform  Generator  (ARB).  This  generator  can  synthesize  an 
approximation  to  any  periodic  waveform.  It  subdivides  the  period  into  N 
time  segments,  where  N  is  selectable  between  1  and  1024.  The  operator 
programs  into  the  instrument's  memory  the  amplitude  of  the  waveform 
(quantized  to  255  levels)  for  each  time  segment,  by  keyboard  entries.  The 
amplitude  scale  factor  and  the  waveform  period  are  also  operator 
selectable.  The  generator  can  generate  a  stepwise  approximation  to  a 
desired  waveform,  or  it  can  be  operated  with  smoothing  between  the  time 
slots  to  generate  a  piecewise  linear  approximation  to  a  desired  waveform. 

The  experimental  setup  is  shown  in  Figure  6,  in  which  the  ARb  is 
used  to  amplitude  modulate  tiie  third  harmonic  carrier  as  the  Reference 
(weighted!  input  to  an  L'tiF  interference  cancellation  system  (ICS).  The  ICS 
Main  (unwe i giited )  input  was  the  third  harmonic  from  a  test  nonlinearity  to 
which  the  amplitude-modulated  fundamental  was  applied.  The  AM  on  the 
fundamental  interference  w.,.s  sinusoidal,  with  the  same  period  as  that 
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FIGURE  6.  EXPERIMENTAL  SETUP  FOR  TIME  DOMAIN  HARMONIC  ENVELOPE  SYNTHESIS 


generated  by  the  ARB.  This  was  accompl  ished  by  using  a  Wavetek  #171 
Synthesized  Signal  Generator  with  an  external  1  MHz  clock  provided  by  the 
10  MHz  clock  oi  the  AKIj  through  a  divide-by-ten  counter.  The  phase  oi  the 
sinusoid  was  adjusted  by  passing  it  through  the  tunable  Krohn-tiite  #310-Ab 
bandpass  Filter. 


The  ARB  was  operated  with  a  period  oi  1.25  khz,  subdivided  into 
N=40  time  slots.  The  ICS  output  was  monitored  on  the  spectrum  analyzer 
used  as  a  fixed-tuned  (at  300  Mhz)  detector  witli  300  kHz  IF  bandwidth.  The 
spectrum  analyzer  then  provided  a  time  domain  display  of  the  ICS  output; 
its  sweep  was  triggered  from  the  ARB.  The  ARB  wavetorin  was  then  manually 
programmed  point-by-point  to  minimize  the  ICS  output  as  seen  on  the 
analyzer.  The  time  slot  modified  at  each  iteration  was  the  one  for  which 
tlie  ICS  output  was  largest.  Smoothing  was  enabled  on  the  ARB. 

Figure  7  shows  the  spectrum  analyzer  display  of  the  ICS  output 
with  no  cancellation  applied,  that  is,  the  ICS  output  with  only  the  output 
of  the  test  linearity  connected  to  the  ICS.  Recall  that  the  analyzer  is 
used  as  a  fixed  tuned  receiver  with  the  horizontal  axis  proportional  to 
time  in  this  experimental  setup. 

Modification  of  the  ARB  waveform  point-by-point  as  described 
above  ai loweo  the  peak  value  of  the  ICS  output  to  be  reduced  by  20  dB. 
Further  .  i  ,  .  t  minimization  could  not  be  accomplished  because  ot  an  ARB 
malfunction  in  which  it  would,  at  random  times,  dump  its  memory.  Further 
experimentation  with  this  equipment  was  aborted  because  of  this 
malfunction,  so  that  more  illustrative  and/or  precise  data  are  not 
avai lable . 

It  was  noted  that  the  ARB  waveform  needed  for  20  db  of 
cancellation  was  not  strictly  symmetrical  about  its  positive  and  negative 
peaks.  This  indicates  that  AM- PM  conversion  effects  were  present,  furtiier 
borne  out  by  the  assymmetry  of  the  spectrum  at  the  ICS  output.  This  is 
very  likely  the  result  of  the  fact  that  other  nonlinearities  were  present 
in  the  experimental  setup  (thus  producing  the  AM-PM  conversion  effects 
described  earlier  in  Section  2.3).  When  the  test  nonlinearity  was  removed, 
the  third  harmonic  interference  decreased  by  about  15  db,  indicating  that 
other  sources  of  distortion  were  present. 

Although  this  laboratory  experiment  did  not  yield  definitive 
data,  it  did  demonstrate  tiiat  harmonic  cancellation  by  synthesis  of  the 
harmonic  envelope  is  a  viable  concept. 

3.3  FOURIER  EXPANSION  HARMONIC  SYNTHESIS 

An  additional  experiment  was  conducted  to  further  evaluate  the 
concept  ot  harmonic  cancellation  by  syntnesis  of  the  harmonic  envelope.  In 
this  experiment  a  sinusoidal  fundamental  envelope  was  used  and  the  harmonic 
envelope  was  synthesized  as  the  sum  of  amplitude  weighted  and  phase  shifted 
harmonics  oi  the  fundamental.  This  harmonic  envelope  synthesis  approach 
relies,  on  the  fact  that  the  harmonic  envelope  is  a  sum  of  powers  of  the 
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fundamental  envelope  as  shown  by  Equation  (2-8).  The  synthesized  harmonic 
envelope  used  in  the  experiment  is 


7 

V  (t)  =  £  /Vcos(27rkft+0  )  (3-1) 

n  k=0  *  * 

while  the  amplitude  modulation  of  the  fundamental  carrier  is  cos(2  ft). 
The  and  0,  in  the  expansion  for  V  (t)  above  are  the  amplitude  and 

phase,  respectively,  of  Lhe  k-th  harmonic  of  the  lundamental  and  were 
controlled  by  an  operator  in  order  to  achieve  the  best  possible 
cancellation . 

Figure  8  shows  the  configuration  of  the  test  setup  used  for  this 
experiment.  The  upper  channel  provides  a  sinusoidally  modulated  100  MHz 
carrier  which  is  driven  into  a  nonlineaity  to  create  a  third  harmonic 
distortion  product.  The  sinusoidal  amplitude  modulation  for  the  100  Miiz 
fundamental  is  generated  by  a  computer  driving  a  D/A  converter.  This  is 
described  in  more  detail  below.  The  third  harmonic  output  of  the 
nonlinearity  is  used  as  the  input  to  the  unweighted  channel  of  the  MX-500 
ICS. 


The  lower  channel  provides  the  synthesized  third  harmonic  by 
first  gcieriting  a  300  MHz  carrier  through  a  bandpass  limiter  and  then 
applying  the  synthesized  third  harmonic  envelope  by  amplitude  modulating 
the  300  MHz  carrier.  The  harmonic  modulation  is  derived  from  the  computer 
as  is  the  fundamental  modulation.  The  resulting  synthesized  third  harmonic 
amplitude  modulated  carrier  is  used  as  an  input  to  the  weighted  channel  of 
the  MX-500  ICS.  The  ICS  is  used  to  effect  the  cancellation  of  the  third 
harmonic  interference  from  the  nonlinearity  by  appropriately  controlling 
the  amplitude  and  phase  of  the  synthesized  third  harmonic.  The  ICS  output 
is  monitored  by  use  of  a  spectrum  analyzer  to  measure  the  cancellation. 
The  waveshape  of  the  synthesized  third  harmonic  amplitude  modulation  is 
varied  by  operator  control  to  achieve  the  greatest  cancellation  at  the  ICS 
output . 


In  order  to  generate  the  fundamental  and  the  synthesized 
amplitude  modulation  waveforms  with  the  absolute  phase  coherence  required 
and  to  permit  the  amplitude  and  phase  of  the  harmonics  in  the  synthesized 
modulation  waveform  to  be  independently  controlled,  these  waveforms  were 
digitally  generated  by  a  computer  and  D/A  converter  for  use  in  the  balanced 
modulators.  The  computer  operated  in  two  different  modes  in  order  to 
generate  these  waveforms:  the  compute  mode  and  the  output  mode.  In  the 
compute  mode  the  computer  accepts  commands  from  the  operator  via  the  CRT 
terminal.  These  commands  established  the  values  of  A^  and  8^  for 
synthesis  of  the  harmonic  modulation  waveform.  After  the  operator  entered 
A^  and  6^  for  all  k  =  0,1,. ..,7  (this  operation  was  simplified  by 

entering  changes  only  in  order  to  expedite  the  optimization  of  the 
synthesized  envelope),  the  computer  generated  two  tables  of  data  each  32 
words  long  and  stored  them  in  memory.  One  table  contained  the  values  for 
the  fundamental  sinusoidal  modulation  (cos27rft)  while  the  second  table 
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FIGURE  8.  EXPERIMENTAL  SETUP  FOR  FOURIER  EXPANSION  HARMONIC  ENVELOPE  SYNTHESIS 


contained  the  values  of  the  synthesized  harmonic  V(t)  as  given  in  (3-1). 


The  second  mode  of  operation  commenced  after  the  above  tables 
were  formed.  In  this  mode  the  computer  repetitively  updated  each  of  the 
two  D/A  converters  with  data  from  the  appropriate  table  in  memory,  thus 
providing  analog  waveforms  at  the  D/A  converter  outputs  corresponding  to 
the  data  stored  in  each  table.  The  tables  were  read  out  repetitively  by 
returning  to  the  first  entry  after  the  last  entry  was  used,  thereby 
producing  periodic  waveforms.  The  fundamental  frequency  and  hence  the 
frequency  of  the  first  seven  harmonics  generated  in  this  way  was 
established  by  control  of  the  frequency  of  the  external  clock  supplied  as 
an  input  to  the  D/ A  converters.  In  order  to  assure  sufficient  precision  in 
the  generation  of  the  waveforms,  system  parameters  and  programs  were 
established  to  provide  at  least  four  samples  per  period  of  the  highest 
frequency  used,  the  seventh  harmonic  of  the  fundamental  modulation.  The 
D/A  converters  were  updated  at  a  rate  of  32f,  where  f  is  the  frequency  of 
the  fundamental  modulation.  With  these  parameters,  the  fundamental 
frequency  could  be  as  high  as  940  Hz.  The  D/A  converters  have  a  resolution 
of  12  bits.  Figure  9  is  a  photograph  of  the  experimental  setup. 

The  fundamental  modulation  frequency  f  has  to  be  selected 
carefully  within  the  following  constraints.  First,  it  must  be  a  frequency 
low  enough  so  that  the  D/A  converter  clock  can  run  at  32f.  Second,  it 
should  be  faster  than  the  tracking  speed  of  the  MX-500  ICS  so  that  the  ICS 
does  not  follow  the  AM.  Third,  it  should  be  low  enough  so  that  mismatches 
of  the  tunable  UHF  BPF's  (approximately  500  kHz  3  dB  bandwidths)  do  not 
limit  the  achievable  cancellation.  Fourth,  f  should  exceed  300  Hz,  so  that 
the  spectral  lines  are  resolvable  by  a  spectrum  analyzer  with  a  minimum 
resolution  of  300  Hz. 

The  procedure  for  manually  selecting  the  values  of  A^  and  0^  is 
based  on  the  acquisition  process  that  an  adaptive  system  would  follow  if 
each  harmonic  were  adaptively  controlled  by  the  LMS  algorithm.  The 
amplitude  corresponding  to  the  strongest  sideband  in  the  error  signal  at 
the  ICS  output  is  the  one  whose  value  should  be  changed  most  during  the 
next  LMS  iteration.  Hence,  the  manual  adjustment  procedure  at  each 
iteration  is  to  vary  the  amplitude  and  phase  corresponding  to  the  strongest 
sideband  in  the  error  signal.  This  applies  to  all  harmonics  except  the 
fundamental,  which  should  be  left  at  its  maximum  value  since  the  MX-500 
will  adaptively  adjust  it.  Therefore,  the  nulling  procedure  involves  a 
trial  and  error  approach  of  modifying  the  harmonic  whose  contribution  to 
the  ICS  error  output  is  greatest,  continually  seeking  the  set  of 
A^,  0^,  k  =  0,1,.. .,7  which  provides  the  smallest  residue. 

Figure  10  shows  the  results  of  an  experiment  using  the  above 
procedure.  The  fundamental  modulation  frequency  is  781  Hz.  The  upper 
photograph  shows  the  spectrum  of  the  ICS  output  with  no  cancellation 
(weighted  channel  disabled).  The  carrier  frequency  is  300  Mhz  and  the 
sideband  structure  caused  by  the  nonlinear  third  harmonic  distortion  is  as 
shown.  The  lower  photograph  shows  the  cancelled  ICS  output  after  the 
adjustment  procedure  described  above  has  been  executed.  The  cancellation 
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SPECTRUM  OF  THE  ICS  OUTPUT 


for  the  central  sidebands  is  30  dB  for  the  component  at  -900  Hz  (with 
respect  to  band  center),  while  the  sideband  at  +900  H.  is  cancelled  by  35 
dB.  The  sidebands  at  +2700  Hz,  +4500  Uz  and  +6300  riz  are  cancelled  to  a 
level  below  the  residue  at  -900  Hz. 

Note  that  sidebands  at  +8100  Hz  and  beyond  are  not  cancelled  at 
all  since  the  synthesized  harmonic  envelope  only  includes  terms  up  to  the 
seventh  harmonic  (6300  Hz).  Observe  also  that  the  spectrum  of  the 
cancelled  residue  is  not  symmetrical  about  the  center  frequency.  This  is 
due  to  AM- to- PM  conversion  effects  in  the  test  setup  causing  the  two 
spectra  at  the  ICS  inputs  to  have  slightly  different  sideband  phases.  The 
AM-to-PM  conversion  is  caused  by  the  fact  that  the  test  setup  contains  more 
than  one  nonlinearity.  This  effect  is  analyzed  in  Section  2.3. 

Figure  11  shows  the  uncancelled  and  cancelled  ICS  output  for  the 
same  experimental  conditions  as  in  Figure  10  except  that  the  analyzer 
resolution  is  3  kHz.  This  display  provides  a  more  integrated  measure  of 
ICS  output  and  corresponds  co  what  would  be  present  in  a  3  kHz  wide 
receiver  tuned  to  some  part  of  the  interference  spectrum.  The  30  dB 
cancellation  at  band  center  is  clearly  shown.  Figure  12  shows  the  waveform 
of  the  synthesized  harmonic  modulation. 

Table  I  is  a  tabulation  of  the  amplitude  and  phases  of  the 
harmonic  components  ot  the  synthesized  harmonic  modulation.  The  0.001 
amplitude  value  is  the  smallest  amplitude  that  could  be  entered.  Since 
this  is  more  than  60  db  below  the  fundamental,  these  components  have  no 
effect  on  the  results. 

Table  1  -  Harmonic  Components  of  Synthesized  Harmonic  Modulation 

Phase  9^  (°) 
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0.001 
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5.495 
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0.001 
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3 

1.757 

-4 

4 

0.001 

0 

5 

0.446 

180 

6 

0.001 

0 

7 

0.1122 

180 

Harmonic  Number  (k)  Amplitude 
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FIGURE  12 

WAVEFORM  OF  SYNTHESIZED  HARMONIC  MODULATION 
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DISTORTION  PRODUCT  CANCELLATION  SYSTEM  CONFIGURATION 


In  this  section  the  configuration  of  systems  designed  for  the 
cancellation  of  harmonic  and  intermodulation  distortion  are  described. 
These  systems  rely  on  some  means  of  synthesizing  the  envelope  of  the 
distortion  products  (DP)  involved  from  a  knowledge  of  the  fundamental 
envelopes.  Two  basic  techniques  for  this  synthesis  are  described.  We 
begin  by  considering  the  cancellation  of  harmonic  interference  and  apply 
the  same  concpets  to  the  more  general  problem  of  intermodulation  products 
(IMP). 

4.1  HARMONIC  ENVELOPE  SYNTHESIS 

Typically,  the  harmonic  interference  coupled  into  a  receiver  is 
due  to  only  one  of  the  harmonics  of  the  offending  transmitter.  In  this 
case  the  interference  is  given  by  (2-8)  for  an  appropriate  selection  of 
harmonic  number  n.  The  objective  of  a  harmonic  interference  cancellation 
system  is  to  cancel  the  interfering  harmonic  signal  to  an  acceptably  low 
level  by  making  use  of  the  transmitted  fundamental  signal.  Figure  13  shows 
the  configuration  of  a  system  which  cancels  harmonic  interference.  The 
system  operates  by  modulating  a  harmonic  carrier  with  a  synthesized 
harmonic  envelope.  The  resulting  synthesized  harmonic  is  applied  to  the 
weighted  channel  of  a  standard  ICS  to  provide  cancellation.  The  ICS 
provides  the  amplitude  and  phase  control  of  the  synthesized  harmonic  RF  (or 
IF)  signal  at  a  rate  which  is  rapid  enough  to  accommodate  any  changes  in 
the  coupling  path  from  the  nonlinearity  to  the  receive  antennas. 

The  harmonic  envelope  synthesizer  operates  on  the  envelope  of  the 
fundamental  signal  to  create  the  envelope  of  the  received  harmonic.  It  is 
directed  by  the  error  signal  from  the  ICS.  The  adaption  time  of  the 
synthesizer  is  slow.  Once  adapted,  it  need  respond  only  to  changes  in  the 
characteristics  of  the  nonlinearity  which  is  assumed  to  be  basically 
nonvarying. 

The  system  of  Figure  13  is  adequate  to  cancel  the  harmonic 
distortion  produced  by  one  nonlinearity  but  will  be  unable  to  provide 
cancellation  for  the  same  harmonic  produced  by  multiple  nonlinearities 
(which  need  not  be  the  same)  excited  by  the  same  fundamental  signal.  As 
shown  in  Section  2.3,  harmonic  energy  receive  from  multiple,  memoryless 
nonlinearities  exhibits  phase  modulation  which  is  dependent  on  the 
amplitude  modulation  of  the  fundamental.  This  AM— to— PM  conversion  requires 
a  more  complex  system  to  effect  cancellation  of  the  harmonic. 

Figure  14  shows  the  configuration  of  a  system  for  multiple 
nonlinearity  generated  harmonic  interference  cancellation.  As  can  be  seen, 
this  system  contains  a  synthesizer  for  both  the  inphase  and  quadrature 
components  of  the  received  interference.  The  inphase  and  quadrature 
components  of  the  ICS  output  are  used  as  feedback  to  the  appropriate 
synthesizer.  This  system  synthesizes  an  inphase  and  quadrature  harmonic 
envelope  and  combines  them  exactly  as  Equation  (2-18)  indicates,  thus 
matching  the  structure  of  the  received  interference. 
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The  crux  of  both  systems  described  above  is  the  synthesis  of  an 
harmonic  envelope  from  knowledge  of  the  fundamental.  The  residue  from  the 
ICS  is  available  as  feedback  to  guide  the  synthesis  procedure.  Means  of 
implementing  this  synthesis  operation  are  described  in  the  following. 

In  the  discussion  which  follows,  x(t)  is  used  to  denor  -  the 
fundamental  envelope.  The  first  step  in  the  synthesis  of  the  envelope  is 
the  expansion  of  the  fundamental  envelope  x(t)  in  a  convenient  set  of  basis 
functions  g.(x).  The  selection  of  the  g.(x)  is  discussed  later.  A 
weighted  sum  of  these  basis  functions  is  used  to  synthesize  the  envelope  of 
the  harmonic  which  we  denote  by  y(t).  Thus  the  harmonic  envelope  is  given 
by 

y (t )  =  yU(t)]  =  l  cigi[x(t)]  (4-1) 


where  the  c  are  constants  and  the  specific  dependence  of  y(t)  on  x(t)  is 
expressed  Dy  the  first  equality.  Note  that  for  the  system  shown  in  Figure 
14,  two  synthesized  envelopes  y-^t)  and  y(.(t)  are  needed.  The  value  of 
the  c  in  the  expansion  for  y(t)  above  must  be  such  as  to  minimize  the 
power  out  of  the  ICS. 

Figure  15  shows  the  operations  required  to  implement  the  above 
processing.  The  fundamental  envelope  is  expanded  in  the  set  of  basis 

functions,  g.(x)  and  each  of  these  expansion  functions  is  then  multiplied 
by  a  constant  c^  and  summed.  The  figure  shows  the  processing  required 
for  envelope  synthesis  in  a  cancellation  system  which  has  a  single  channel 
as  shown  in  Figure  13.  The  summation  of  t he  weighted  basis  functions  which 
is  the  synthesized  envelope  is  used  to  modulate  the  harmonic  carrier. 

Adaptive  control  of  the  c  is  implemented  by  an  LMS  feedback 
loop  which  acts  to  minimize  the  output  of  the  canceller.  The  LMS  control 
loop  is  well  suited  to  this  application  since  it  acts  to  minimize  the 
output  power  of  the  canceller.  The  ICS  output  is  first  correlated 

(synchronously  detected)  with  the  synthesized  carrier  to  extract  the  error 
at  baseband  and  then  is  correlated  with  g.(x)  to  determine  the  error 
component  associated  with  g^(x).  This  error  signal  is  applied  to  the 
multiplier  through  a  lowpass  filter. 

The  means  of  implementing  the  various  signal  processing 

operations  shown  on  Figure  15  depends  on  the  specific  form  selected  for  the 

g.(x).  Several  techniques  of  adaptively  synthesizing  the  harmonic 
envelope  have  been  identified.  The  approaches  to  envelope  synthesis  fall 
in  two  categories:  gate  functions  and  polynomial  expansions. 

4.1.1  Gate  Function  Synthesis 

The  use  of  gate  functions  in  envelope  synthesis  is  a 

straightforward  approach  well  suited  to  digital  processing  circuits.  Gate 
functions  are  binary  valued  functions  which  are  useful  for  defining  when 
variable  is  within  certain  bounds.  Formally,  a  gate  function  g.(x)  is 


a 


defined  as  follows: 


g^(x)  =  1,  if  x1_1  <  x  -  xi 
=  0,  otherwise 


(4-2) 


for  i  =  1,2 .  Figure  16  shows  several  gate  functions.  They  can  be  used 

to  expand  an  arbitrary  function  f(x)  as  follows: 

N 

f(x)  =  I  c  g,(x)  (4-3) 

k=l  k  * 

which  is  a  stepwise  function  of  x  which  can  assume  an  arbitrary  shape 
depending  on  the  selection  of  c^.  For  any  x  in  an  allowed  range,  one  and 
only  one  g^  has  nonzero  value  and  f(x)  assumes  the  value  c^  for  that  x. 

For  our  application,  x  represents  the  fundamental  envelope  and  f(x) 
represents  the  harmonic  envelope.  A  typical  f(x)  is  shown  in  Figure  17. 
The  objective  of  adaptive  harmonic  envelope  synthesis  is  to  select  the  set 
of  c^  so  that  the  system  output  is  minimized. 

A  method  is  given  in  ill]  for  computing  each  c  coefficent  to 
minimize  mean-square  error  e7"  between  the  approximated  output  and  a  desired 
output  z(t).  The  value  of  c  should  be  set  so  that  9e^79c  =  0.  This  is 

achieved  when  the  error  n 

N 

e  =  z(t)  -  J  cngn(x)  (4-4) 

n=l 


is  orthogonal  to  each  gate  function,  g^(x).  Thus, 
Egn00  =  [z(t)-cngn(x)]gn(x) 

=  [z(t)-cn]gn(x) 


where  the  overbar  represents  the  averaging  process  used  in  defining  the 
mean  square  error.  If  a  time  average  is  used,  then  adjustment  of  the 
coefficients  {c^}  to  satisfy  (4-5)  in  steady  state  may  be  implemented  by 
the  LMS  (Least  Mean  Squares)  feedback  algorithm.  Each  coefficient  is 
adjusted  by  a  high  gain  negative  feedback  loop,  which,  in  the  steady  state, 
satisfies  Equation  (4-5).  The  LMS  algorithm  generates  an  incremental 
negative  feedback  correction  to  each  c that  is  proportional  to  the 

[11]  M.  Schetzen,  "Determination  of  Optimum  Nonlinear  Systems  for 
Generalized  Error  Criteria  Based  on  the  Use  of  Gate  Functions,"  IEEE  Trans 
on  Information  Theory,  January  1965,  pp.  117-125. 
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derivative  of  e  with  respect  to  c^,  that  is: 


dc  (t)  a  2 

TT-  •  t-  '  -2u£®n(x) 

n 


Thus,  in  a  sampled-date  control  system 

c  (j+1)  =  c  (j)  -  2ye(j)g  [x(j>] » 
n  n  n 

while  in  a  continuous  system 


c 

n 


(t)  =  -2y  /  e(T)g  fx(T)] 
0 


dt 


(4-6) 


(4-7) 


(4-8) 


Equation  (4-7)  is  particularly  easy  to  implement  digitally.  It 
tells  us  that  whenever  the  n-th  gating  function  is  activated,  then  c  is 
incremented  by  -2y e  from  its  previous  value.  When  the  n-th  gating  function 
is  not  actuated,  then  c  remains  at  its  previous  value.  Since  the  gating 
function  boundaries  (ft  >xn+j)  do  not  overlap,  only  one  gating 
function  is  activated  at  any  one  time.  Thus,  the  weighted  sum  of  gating 
functions  used  in  (4-4)  to  approximate  z(t)  reduces  to 

z(t)  =  cr  when  xR<x(t)  -  xn+J  and  e-K)  (4-9) 

If  we  apply  x(t)  to  an  analog-to-digital  converter  that  generates  the 
number  n  when  x  <x(t)<xn+j,  then  that  number  may  be  used  to  address  a 
memory  in  which  tRe  c  ' s  are  stored.  The  addressed  c  is  used  to 

approximate  z(t)  ana  is  also  updated  to  a  new  value  byn incrementing  it  by 
-2ye. 


Figure  18  depicts  the  processing  required  to  implement  harmonic 
envelope  synthesis  using  gate  functions.  The  analog-to-digital  (A/D) 
converted  fundamental  envelope  is  used  as  an  address  for  the  random  access 
memory.  This  selection  of  a  unique  address  for  each  unique  code  out  of  the 
A/D  converter  (corresponding  to  a  range  of  envelope  values)  implements  the 
gate  function  expansion  described  above.  The  number  stored  in  the  address 
represents  the  value  c^  assigned  to  g^(x).  The  contents  of  the  current 
address,  c^,  are  read  out  of  memory  and  digital-to-analog  (D/A)  converted 
to  form  the  input  to  the  modulator.  Adaptive  control  of  the  value  of  c, 
is  achieved  by  determining  the  error  in  the  system  output  (by  correlation) 
and  incrementing  the  contents  of  the  current  address  by  a  negatively  scaled 
version  of  this  error.  (This  processing  implements  the  operations  of  (4-7) 
above  which  indicates  that  when  the  c.  reach  steady  state  the  error  is 
small.)  1 
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RECEIVED 


FIGURE  18.  HARMONIC  SYNTHESIS  USING  GATE  FUNCTIONS 


.«  . 


This  implementation  is  seen  to  involve  straightforward  digital 
processing  techniques  where  the  only  arithmetic  operation  involved  is 
addition. 


4.1.2  Polynomial  Expansion  Envelope  Synthesis 

An  alternate  means  of  expanding  the  nonlinearity  to  synthesize 
the  harmonic  envelope  is  to  make  use  of  standard  polynomial  expansions.  A 
wide  variety  of  possibilities  is  available.  The  simplest  is  a  power  series 
expansion  in  which  the  expansion  functions  are  given  by 

g. (x)  =  xi+n,  i  =  0,1,...  ,N 
where  n  is  the  order  of  the  harmonic. 


Other  polynomials  such  as  Tchebychev  are  available  for  expanding 
the  nonlinearity.  In  this  case  the  are  given  by 


8i  = 


where  n  is  the  order  of  the  harmonic  and  T  (x)  is  the  Tchebychev 
polynomial  of  order  i.  Several  of  these  are  listed  below: 


by: 


T0(x)  =  1 
Tj(x)  =  x 
T2(x)  =  2x2-1 
T^(x)  =  4x^-3x 
etc . 

In  either  case  the  expansion  of  the  harmonic  envelope  is  given 
N 

f(x)  =  l  c  g  (x)  (4-10) 

i=0 


where  x  is  the  fundamental  envelope.  As  before,  the  coefficients  c  must 
be  found  to  minimize  the  system  output.  Figure  19  shows  the  system 
configuration  for  utilization  of  polynomials  in  the  synthesis  of  the 
harmonic.  The  envelope  of  the  fundamental  is  A/D  converted  and  the  g  (x) 
computed  for  the  envelope  value  x.  Then  Ec  g  (x)  is  formed  and  used  !o 
modulate  the  harmonic  carrier.  Adaptive  control  of  each  c  is  obtained 
by  correlation  with  the  harmonic  carrier  and  then  correlating  each  g 
with  this  error  signal.  Correlating  each  g  with  the  error  produces 
correction  signals  which  are  applied  to  the  c  to  minimize  the  system 
output.  The  corrections  are  negatively  Icaled  and  lowpass  filtered  to 
implement  the  LMS  control  algorithm. 
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FIGURE  19.  HARMONIC  SYNTHESIS  USING  POLYNOMIALS 


The  response  time  of  envelope  synthesis  can  be  improved  by 
selecting  expansion  polynomials  g.(x)  that  are  orthogonal.  Using 
orthogonal  polynomials  reduces  the  cross  coupling  in  the  correlation 
processing  (ideally  to  2ero)  and  permits  the  coefficients  c  to  be  found 
more  rapidly.  The  Tchebychev  polynomials  are  orthogonal  ana  provide  the 
additional  advantage  that  minimization  of  the  average  power  out  of  the 
canceller  also  minimizes  the  peak  power. 

As  can  be  seen  from  Figure  19,  considerably  more  arithmetic 
processing  is  required  to  synthesize  the  harmonic  envelope  using  polynomial 
expansions  in  comparison  to  the  gate  function  approach  (refer  to  Figure 
18).  This  processing  is  within  the  capabilities  of  today's  digital 
processing  components  and  can  be  done  with  a  combination  of  table  look-up 
and  special  purpose  arithmetic  processors. 

4.2  INTERMODULATION  PRODUCT  SYNTHESIS 

Intermodulation  product  (IMP)  interference  is  produced  when  the 
fundamental  signal  from  two  distinct  transmitters  excite  a  nonlinearity  and 
one  of  the  resulting  distortion  components  is  coupled  to  the  receive 
antenna  of  a  receiving  system.  The  objective  of  an  IMP  cancellation  system 
is  to  cancel  the  interfering  IMP  to  an  acceptably  low  level  by  making  use 
of  only  the  fundamental  signals.  The  structure  of  the  IMP  at  a  particular 
frequency  fR  =  mjfj+m^  where  f.  and  f2  are  the  fundamental 
carrier  frequencies  is  derived  in  Section  2.2.  The  IMP  given  by  (2-15) 
where  the  amplitude  modulation  of  the  IMP  y(t)  is  expressed  as  a  weighted 
summation  of  powers  of  the  fundamental  modulations. 

Figure  20  shows  the  configuration  of  a  system  which  operates  to 
cancel  IMP  interference.  The  structure  of  this  system  is  similar  to  that 
of  the  harmonic  cancellation  system.  In  the  present  case  two  fundamental 
signals  must  be  processed  to  generate  a  carrier  signal  at  the  IMP  frequency 
and  to  synthesize  the  IMP  envelope.  Generation  of  the  IMP  carrier  can  be 
accomplished  by  application  of  standard  hardlimiting  and  filtering 
techniques.  Synthesis  of  the  IMP  envelope,  as  before,  is  the  difficult 
part  of  the  system  design. 

Figure  21  shows  the  structure  of  the  processing  required  for  IMP 
envelope  synthesis  by  a  direct  application  of  Equation  (2-15)  and  its 
associated  definitions  (the  equations  following  (2-15)).  As  these 
equations  show  the  IMP  envelope  is  a  weighted  sum  of  powers  of  the 
fundamental  envelopes  A  (t)  and  A,(t).  The  weighting  coefficients 
C^l»^2  are  related  to  the  coefficients  a  of  the  power  series  expansion 
of  the  nonlinearity.  These  coefficients  must  be  selected  to  closely 
approximate  the  received  IMP  envelope. 

4.2.1  Synthesis  of  the  IMP  Envelope  by  Use  of  Gate  Functions 

The  amplitude  of  the  IMP  at  f  can  be  synthesized  using  the 
gate  functions  described  earlier  in  a  system  configured  as  shown  in  Figure 
22.  The  system  configuration  is  similar  to  that  used  for  harmonic  envelope 
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FIGURE  20.  IMP  CANCELLATION  SYSTEM 


FIGURE  21.  SYNTEHSIS  OF  IMP  AMPLITUDE  MODULATION 


synthesis  (refer  to  Figure  18).  The  output  of  each  fundamental  transmitter 
is  A/D  coverted  and  used  to  address  into  a  random  access  memory.  When  used 
for  harmonic  envelope  synthesis,  the  output  of  a  single  A/D  was  used  to 
address  into  a  particular  memory  location,  thus  implementing  the  gate 
function  processing.  For  the  harmonic  case,  the  memory  was 
one-dimensional,  that  is,  each  A/D  converter  output  selected  one  memory 
location.  For  use  in  the  synthesis  of  IMP  envelopes,  the  memory  is 
two-dimensional,  that  is,  the  A/D  output  of  one  fundamental  envelope 
determines  the  value  of  one  coordinate  and  the  A/D  output  of  the  second 
fundamental  determines  the  value  of  the  second  coordinate.  A  particular 
memory  location  is  determined  by  both  A/D  outputs,  thus  implementing  a 
two-dimensional  gate  function. 

The  contents  of  the  memory  location  addressed  in  this  way  contain 
a  constant  which  is  the  amplitude  of  the  IMP  corresonding  to  the  two 
fundamental  envelopes.  The  value  of  this  constant  is  continually  updated 
by  the  operation  of  an  LMS  control  loop  as  was  done  in  the  harmonic 
envelope  case  and  defined  by  Equation  (4-7). 

Because  of  the  need  to  drive  the  IMP  synthesis  processing  from 
two  fundamental  envelopes,  the  size  of  the  memory  required  will  increase 
substantially.  In  addition,  due  to  the  increased  number  of  memory 
locations,  the  time  required  for  the  synthesis  processing  to  converge  to 
its  steady  state  will  be  greatly  increased  over  the  adaption  time  for  the 
harmonic  envelope  synthesis  case. 

4.2.2  Polynomial  Expansion  IMP  Envelope  Synthesis 

As  was  done  in  the  case  of  harmonic  envelope  synthesis, 
polynomial  expansion  can  be  used  to  synthesize  the  IMP  envelope.  To  see 
how  this  is  done,  we  begin  by  expanding  the  fundamental  envelope  terms  of 
the  IMP  y(t)  derived  earlier  (see  Equation  (2-15))  as  follows: 


=  l  3r(^1)gr(A1) 

42  =  *  WMV 


where  the  coefficients  3^2^)  and  3g(22)  Can  be  computed  exactly  once  the 
form  of  g{(x)  is  .selected.  Using  these  representations  in  y(t)  as  given 

and  interchanging 


in  (2-15),  the  IMP  envelope  derived  in  Section  2.2, 
orders  of  summation,  y(t)  can  be  expressed  as 


m.  m„ 

y(t)  -  Ax  A2  ll  br>s  gr(Ax)  gs(A2) 


(4-11) 


r  s 
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where 


b 

r,s 


Hl=r  S-2=s 


A, 


ww 


Figure  23  shows  the  system  configuration  which  implements  the 
abov<?  processing.  its  structure  is  similar  to  that  used  for  polynomial 
expansion  of  harmonic  envelopes  shown  in  Figure  19.  As  shown  in  Figure  23, 
control  of  the  coefficients  b  is  implemented  by  use  of  an  U1S  control 
loop.  'Hie  correlation  anti  updating  processing  implements  the  control 
algorithm  shown  in  Figure  14. 


The  arithmetic  processing  required  to  implement  the  operations 
shown  in  Figure  23  is  much  greater  than  for  the  gate  function  synthesis 
approach  shown  in  Figure  22.  However,  the  memory  requirements  are  less  and 
the  acquisition  time  should  be  more  rapid.  More  study  is  required  to 
compare  the  performance  and  complexity  of  the  two  IMF  envelope  synthesis 
approaches. 


4.2.3  System  Configuration  for  Multi  p.le  Non  L  inoar  i_t  ies 

As  was  the  case  for  harmonics,  multiple  non  1  inear i t ies  will 
produce  AM-to-FM  conversion  effects  for  IMF.  By  providing  separate 
synthesis  capability  for  1  and  Q  channels,  it  should  be  possible  to 
generate  a  replica  of  the  composite  received  LMF.  More  study  is  required 
to  evaluate  this  approach  and  to  investigate  means  of  reducing  its 
complexity.  Sucli  a  system  is  roughly  twice  the  complexity  of  the  systems 
shown  in  Figure  22  or  23. 
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5.U  SIMULATION  OF  ENVELOPE  SYNTHESIS  BY  GATE  FUNCTIONS 

In  order  to  determine  the  effectiveness  of  an  harmonic 
cancellation  system  which  employs  gate  Function  synthesis  of  the  harmonic 
envelope,  system  operation  was  simulated  on  a  digital  computer  and  the 
performance  of  the  system  evaluated  for  a  variety  of  external 
nonlinearities  and  fundamental  amplitude  modulation  waveforms.  The 
simulation  was  configured  in  such  a  way  as  to  accurately  represent  the 
digital  processing  that  would  be  implemented  in  a  practical  system 
realization.  In  the  following,  the  simulation  is  described  in  detail  and 
the  results  discussed. 

5.1  DESCRIPTION  OF  SIMULATION 

The  harmonic  cancellation  system  simulated  is  that  shown  in 
Figure  14  which  includes  separate  I  and  Q  channels  to  properly  synthesize 
the  harmonic  envelope  for  the  situation  where  there  is  AM-to-PM  coversion 
due  to  the  presence  of  multiple  nonlinearities.  The  I  and  Q  channel 
harmonic  envelope  synthesis  is  implemented  by  use  of  gate  functions  as 
shown  on  Figure  18. 

The  simulation  operates  with  the  complex  envelopes  of  all  RF 
signals  actually  present  in  the  system.  Therefore,  the  ICS  which  produces 
the  cancellation  by  weighting  the  synthesized  harmonic  and  combining  it 
with  the  received  interference  is  not  included  in  the  simulation  of  the 
system.  This  exclusion  is  not  significant  in  terms  of  the  evaluation  of 
the  harmonic  cancellation  system  performance,  since  the  ICS  performance  is 
not  the  limiting  factor  in  this  system.  It  functions  to  adjust  the 
amplitude  and  phase  of  the  synthesized  harmonic  in  order  to  compensate  for 
small  differences  in  the  interference  coupling  path. 

Figure  24  shows  the  structure  of  the  simulation.  The  simulation 
program  was  written  in  FORTRAN  and  run  on  a  PDP-11/35  computer.  Each 
iteration  of  the  program  begins  with  the  generation  of  a  sample  of  the 
fundamental  modulation  waveform.  This  sample  is  used  to  generate  a  sample 
of  the  output  of  at  least  one  nonlinearity.  Up  to  three  nonlinearities 
were  simulated  at  one  time.  All  nonlinearities  were  modelled  as  polynomial 
functions  of  the  fundamental  waveform. 

The  fundamental  sample  is  also  used  to  generate  the  phase  of  the 
synthesized  harmonic  carrier.  The  carrier  phase  is  the  angle  of  the 
(possibly)  complex-valued  fundamental  amplitude  modulation  waveform.  In 
order  to  assure  that  carrier  is  present  for  both  I  and  Q  channels,  a  fixed 
45  phase  offset  was  added  to  the  phase  found  as  described  above.  The 
envelope  detector  function  is  accomplished  by  forming  the  magnitude  of  the 
fundamental.  The  function  of  A/D  conversion  is  implemented,  by  converting 
from  floating  point  to  integer  representation  in  the  simulation  program. 
Prior  to  integer  conversion,  the  envelope  is  clipped  to  represent  the 
effect  of  an  A/D  converter  with  a  limited  number  of  bits. 


GENERATE 
SAMPLE  OF 


FIGURE  24.  HARMONIC  ENVELOPE  SYNTHESIS  SIMULATION 


The  clipper  function  is  the  following,  where  Y  is  the  clipper 

output . 

Y  =  iV0l  if  |VQ|  $  2N-1 

=  2N-1  if  j VQ1  >  2N-1 

This  clipping  function,  when  followed  by  integer  conversion,  represents  the 
operation  of  an  A/D  converter  with  N-bit  output  resolution. 

The  integer  output  of  the  A/D  converter  is  used  as  an  address  to 
select  one  entry  from  each  of  two  tables.  The  two  tables  represent  the  I 
and  Q  components  of  the  synthesized  harmonic  envelope  and  the  addressing 
described  above  implements  the  mathematical  structure  of  the  gating 
function.  The  bit  length  of  each  entry  in  the  tables  was  restricted  to  N 
bits  to  accurately  represent  a  practical  system  implemented  with  fixed 
point  arithmetic  and  N-bit  processing  and  storage  capabilities. 

The  addressed  1  and  Q  table  values  are  D/A  converted  to  floating 
point  representation  and  multiplied  by  the  complex  the  complex  envelope  of 
the  synthesized  harmonic  carrier  (  unity  magnitude)  to  form  the  synthesized 
modulated  harmonic.  The  complex-valued  harmonic  interferencee  is 
subtracted  from  the  synthesized  harmonic  to  form  the  canceller  error 
signal.  After  multiplication  by  the  synthesized  harmonic  carrier  to  effect 
correlation  or  synchronous  detection  and  negative  amplification,  the  error 
is  clipped  and  A/D-converted  to  an  N-bit  integer  number  as  described 
earlier.  The  current  1  and  Q  values  are  added  to  the  resulting  error 
signal  and  the  results  loaded  into  the  tables  as  the  new  I  and  Q  values. 
The  new  I  and  Q  values  are  clipped  at  a  level  of  2-1  to  represent 
storage  limitations  of  N-bit  words. 

The  processing  described  above  implements  the  LMS  control 
algorithm  for  both  the  I  and  Q  values  of  the  synthesized  harmonic  envelope. 
When  a  particular  entry  in  the  I  and  Q  tables  is  addressed,  the  updating 
equation  for  those  values  is 

X(k+1)  =  X(k)-Ge 

where  X(k+1)  =  new  value 
X(k)  =  old  value 

e  =  error  at  the  canceller  output 
G  =  gain  in  error  path 

The  above  disregards  the  operation  of  the  clipper  and  A/D  converter. 

The  simulation  program  included  the  definition  of  the  fundamental 
waveform  and  the  number  and  form  of  the  nonlinearities.  At  run  time  the 
following  parameters  could  be  selected: 


Number  of  bits  (N) 

Sample  rate 
Gain  (G) 

Fundamental  waveform  frequency 
Run  time 

The  output  data  from  the  program  is  the  time  history  of  cancellation 
defined  as 

C  -  10  log(|e|2/jvJ2) 

where  |e | 7  *  average  of  the  squared  error  at  canceller  output 

| vr|  2  **  average  of  the  squared  harmonic  interference 

The  value  of  | v  | 2  is  also  provided  as  an  output.  The  average  is  computed 
over  one  hundre9  samples  of  the  respective  signals. 

An  additional  output  is  available  on  an  optional  basis.  A  plot 
of  the  values  of  the  I  and  Q  synthesized  harmonic  tables  can  be  obtained. 

Generally,  the  maximum  levels  of  both  the  fundamental  waveform 
and  the  harmonic  interference  are  set  to  utilize  nearly  the  full  dynamic 
range  provided  by  the  number  of  bits  N.  That  is.  the  maximum  values  were 
selected  to  be  approximately  0.9  or  0.95  times  2  . 

Details  of  the  simulation  are  available  from  the  FORTRAN  program 
listing  in  Appendix  B. 

5.2  SIMULATION  RESULTS 

The  simulation  program  was  run  for  a  variety  of  fundamental 
waveforms,  nonlinearities,  numbers  of  nonlinearities  and  canceller  system 
parameters.  The  results  of  these  various  runs  indicate  that  the  gate 
function  synthesis  approach  is  capable  of  creating  the  harmonic  envelope  by 
making  use  of  knowledge  of  the  fundamental  amplitude  modulation. 

A  variety  of  fundamental  amplitude  modulations  were  used.  The 
following  are  typical  examples. 

Amplitude  modulation  (AM):  A(l+asincot) ,  a<l 
AM/Suppressed  carrier  (SC):  Asinyt 
Multi-tone  (MT)  AM/SC:  AjSin^jt  +  A2sinu)2t 

Generally  the  fundamental  modulation  was  selected  to  have  a  maximum 
amplitude  less  than  unity.  The  result  was  then  scaled  prior  to  entry  into 
the  envelope  detector.  The  scaling  constant  is  2-1 ,  the  dynamic  range 
of  the  A/D  converter  for  any  value  of  N. 
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The  same  scaling  is  applied  to  the  output  of  the  composite 
nonlinearity  with  the  nonlinearity  characteristic  selected  to  yield  a 
maximum  envelope  ot  less  than  unity.  This  scaling  assures  that  the  full 
dynamic  range  of  the  N-bit  words  in  the  1  and  Q  envelope  synthesis  tables 
is  utilized  without  overflow.  The  following  nonlinearities  were  used  in 
various  combinations  for  simulation  runs. 


A. 


Diode  Approximation 


V3  = 


u  *  i-t i  +  feu  + 


where  v  *  fundamental  amplitude 

Vj  =  third  harmonic  amplitude 
This  approximation  is  from  reference  [6]. 

3 

B.  Third-Order  Polynomial:  v^  =  A^v  +  A^v 

C.  Fifth  Power  of  Fundamental:  v^  =  Av"* 

These  nonlinearities  represent  distributed  nonlinearities  as 
described  in  Section  2.3.  Therefore,  the  harmonic  amplitudes  of  each 
separate  nonlinearity  are  combined  with  different  RF  phases  to  produce  a 
complex,  composite  harmonic  interference  at  the  input  to  the  canceller. 
All  simulations  were  formally  run  for  the  third  harmonic  case  although  the 
results  are  general  and  apply  to  higher  order  harmonics. 


Table  II  contains  a  summary  of  the  results  of  typical  simulation 
runs  for  a  variety  of  simulation  parameters.  The  fundamental  modulation 
and  nonlinearities  have  been  defined  above. 


Multiple  entries  in  the  nonlinearity  column  indicate  that  those 
nonlinearities  were  present  simultaneously  to  produce  a  composite 
nonlinearity.  The  column  entitled  "Special  Conditions”  indicates  the 
deviations  from  the  standard  conditions  made  during  certain  of  the 
simulation  runs.  Average  harmonic  power  is  the  mean  square  value  of  the 
harmonic  interference  at  the  canceller  input  with  a  0  dB  reference  of  1 
watt.  Average  steady  state  (SS)  cancellation  is  given  by 


C  =  10  log(je|2/jv|2) 


where  | y | 2  =  mean  square  harmonic  interference 
1ST-  mean  square  canceller  output 

The  final  column,  the  acquisition  time,  is  a  measure  of  the  time 
taken  for  the  cancellation  to  reach  a  level  within  2  to  A  dB  of  the  steady 
state  cancellation.  The  values  listed  in  this  column  are  not  precise  and 
should,  be  taken  as  an  approximate  measure  of  the  time  for  the  harmonic 
envelope  synthesis  to  settle. 
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Study  of  simulation  output  data  indicates  that  cancellation 
starts  at  0  dB  and  improves,  finally  reaching  the  indicated  steady  state 
value.  However,  the  improvement  is  not  smooth  and  continuous.  The 
cancellation  is  observed  to  increase  steadily  for  a  time  period  and  then 
decrease  for  a  short  interval,  returning  to  better  cancellation  after  the 
sharp  decrease. 

The  reason  for  this  is  the  fact  that  a  particular  entry  in  the  I 
and  Q  synthesis  tables  is  updated  only  at  those  times  when  the  fundamental 
envelope  falls  in  the  range  of  signal  levels  which  yield  the  proper 

address.  During  the  adaption  process  certain  table  entries  are  skipped 

over  more  than  others  and,  therefore,  remain  in  their  initial  zero  or 
partially  adapted  state  for  a  longer  time.  The  time  at  which  they  begin  to 
adapt  or  resume  adaption  occurs  when  the  fundamental  is  sampled  at  the 
correct  time  to  produce  the  required  address.  This  process  can  lead  to 
situations  where  the  adaption  of  certain  table  entries  takes  considerably 
longer  than  other  entries,  thus  producing  bursts  of  poor  cancellation 
during  the  adaption  process. 

Another  factor  affecting  the  cancellation  data  is  the  way  in 

which  they  are  computed.  The  mean  square  canceller  output  is  found  by 

arithmetically  averaging  the  squared  magnitude  of  the  error  signal  over 
blocks  of  100  samples.  This  averaging  over  blocks  also  tends  to  produce 
somewhat  nonuniform  results. 

The  acquisition  time  listed  in  Table  II  is  the  approximate  time 

it  takes  for  the  cancellation  to  remain  within  2-4  dB  of  the  final  value 

with  no  short  intervals  of  poor  cancellation. 

As  the  data  in  Table  II  indicate,  the  harmonic  envelope  synthesis 
can  be  accomplished  very  accurately.  For  those  cases  where  the  gain  is 
unity,  the  steady  state  cancellation  is  about  equal  to  the  average  harmonic 
power.  The  limitation  to  ultimate  cancellation  in  these  cases  is  the 
limited  number  of  I  and  Q  table  entries  (2  -1  locations)  and  the 
limitation  of  each  of  these  to  N  bits.  These  two  effects  appear  as  a 
minimum  resolution,  placing  a  lower  bound  on  the  achievable  cancellation. 
Note,  that  the  case  where  G  “  0.25  (entry  7)  achieves  considerably  reduced 
cancellation. 

The  acquisition  time  for  those  cases  run  under  standard 
conditions  is  between  1  and  2  seconds.  For  entries  4  and  5,  the 
acquisition  time  is  greatly  reduced  from  this,  because  there  are  many  fewer 
I  and  Q  synthesis  table  entries  to  update.  The  acquisition  time  for  entry 
7  where  the  gain  is  reduced  is  smaller  than  those  cases  run  with  standard 

conditions  because  the  steady  state  cancellation  is  lower.  The  1.2  second 

acquisition  time  for  entry  9  which  is  for  the  multi-tone  AM/SC  case 
indicates  that  variations  in  fundamental  envelope  waveform  affect  the 
settling  time  of  the  synthesis  procedure. 

The  values  of  the  I  and  Q  synthesis  tables  are  plotted  in  Figure 
25  for  the  simulation  run  of  entry  6  in  Table  II.  For  this  simulation  run, 
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the  RF  phases  of  the  A  and  B  nonlinearities  were  selected  to  align  exactly 
with  the  phases  of  the  I  and  Q  components  of  the  synthesized  harmonic 
carrier,  respectively.  The  result  of  this  alignment  is  that  the  I  channel 
adapts  to  synthesize  nonlinearity  A  while  the  Q  channel  adapts  to 
synthesize  the  B  nonlinearity.  Figure  25  is  a  plot  of  the  normalized  value 
of  A  and  B  synthesis  tables  as  a  function  of  the  address  of  the  tables. 
The  plot  is  normalized  to  the  maximum  value  in  either  the  I  or  Q  table. 
Each  curve  represents  the  nonlinear  transfer  function  of  the  appropriate 
nonlinearity. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


b  .0 


This  study  was  directed  at  the  analysis  of  the  characteristics  of 
harmonic  interference  and  the  synthesis  and  evaluation  of  techiques  to 
effect  cancellation  of  such  interference.  The  results  were  extended  to 
indicate  approaches  to  cancellation  of  more  complex  distortion  product  (DP) 
interferences  including  harmonics  from  multiple  nonlinearities  and 
intermodulation  products  caused  by  two  fundamental  signals.  During  the 
initial  stages  of  the  investigation,  a  survey  of  the  literature  was 
conducted  in  order  to  determine  if  any  existing  work  was  available  to 
assist  in  the  problem  formulation  or  solution.  In  spite  of  the  current 
renewed  interest  in  nonlinear  system  analysis,  no  work  was  uncovered  which 
addressed  the  problem  from  the  standpoint  of  detailed  harmonic  or  DP 
waveform  analysis;  hence  the  available  work  is  not  directly  applicable  to 
the  DP  cancellation  problem. 

In  addition  to  the  analysis  of  harmonic  waveforms,  extensive 
laboratory  experiments  were  conducted  to  evaluate  the  extent  to  which 
typical  nonlinearities  are  memoryless.  The  conclusion  of  these  experiments 
is  that,  at  least  to  a  level  corresponding  to  20  to  30  dB  of  cancellation 
over  limited  bands,  the  nonlinearities  are  memoryless.  These  experiments 
also  provided  initial  verification  of  the  concept  of  harmonic  envelope 
synthesis  for  harmonic  interference  cancellation. 

Further  evaluation  of  harmonic  envelope  synthesis  by  use  of  gate 
function  processing  was  obtained  by  simulating  the  essential  parts  of  a 
harmonic  cancellation  system.  The  results  of  the  simulation  indicate  that 
the  approach  is  viable  and  realizable  with  available  analog  and  digital 
signal  processing  components.  The  complexity  of  such  an  implementation 
appears  to  be  reasonable. 

The  results  of  the  study  indicate  that  an  harmonic  cancellation 
system  employing  the  concept  of  harmonic  envelope  synthesis  is  feasible  and 
that  the  next  step  should  be  the  design  and  implementation  of  an 
experimental  canceller.  Many  practical  design  issues  remain  to  be  resolved 
prior  to  implementation.  Among  them  are  the  following. 

The  effect  of  variations  in  the  coupling  between  the 
fundamental  transmitter  and  the  harmonic  producing 
nonlinearities.  Such  variations  may  require  more  rapid 
synthesis  adaption  or  greater  dynamic  range. 

Coupling  variations  between  the  nonlinearity  output  and  the 
victim  receiver  may  be  completely  accommodated  by  the  ICS  or 
may  need  additional  level  control  processing. 

The  effect  and  relationships  between  various  fundamental 
waveforms  and  the  synthesis  update  rate. 

A  consideration  of  the  total  RF  bandwidth  that  must  be 
processed  at  the  victim  receiver  input  to  properly  represent 
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Optimization  of  the  number  of  bits  and  the  memory  size 
requirements  for  a  desired  level  of  cancellation.  This 
should  include  a  consideration  of  nonlinear  A/D  and  D/A 
conversion  characteristics  (companding)  to  make  most 
effective  use  of  the  available  dynamic  range. 

With  these  system-related  questions  resolved,  the  design  of  an  experimental 
model  can  proceed.  Of  course,  many  design  issues  remain  to  be  resolved, 
but  these  are  the  typical  implementation  questions  associated  with  any 
canceller  system  design. 


APPENDIX  A 


CANCELLATION  ACHIEVABLE  BY  DIFFERENT  NONLINEARITIES 
A1.0  INTRODUCTION 

As  derived  in  Section  2,  the  waveform  at  the  output  of  a  nonline¬ 
arity  characterized  by  the  power  series  expansion 


VOUT(t)  =  l  VlN(t) 
k=0 


is  given  by 


OUT 


00  a  T(k+l)V  (t)cos[nL0  t+n(j)(t:)] 

(t)  =  i  -  c 

k=n 


2kr(^4i)r(~+i) 


(A-l) 


(A-2) 


(k-n  even) 

where  vTxl(t)  =  V(t)  cos  w  t+cf>(t)  .  (A-3) 

IN  c 

Thus  the  principal  envelope  term  at  the  n-th  harminic  is  (k=n  in  sum) 


e(t)  = 


anVnu) 


(A-A) 


and  it  is  corrupted  by  the  higher  order  terms  (k  =  n+2,n+3,...) 

+  (n+6)(n+5)(n+A)  a  vn+6(t)  +  ...  (A-5) 

(2n+?) (3)  n+6 

Assuming  that  the  corrupting  terms  are  uniformly  convergent,  the  principal 
corrupting  term  is 


,  .  (n+2)  n+2,  . 

El(t)  "  2n+2  an+2V 


(A-6) 


63 


A2.0  ACHIEVABLE  CANCELLATION  IN  TERMS  OF  NONLINEARITY  EXPANSION 

COEFFICIENTS 

Let  us  model  the  cancellation  problem  as  shown  in  Figure  A-l.  We 
want  to  find  the  achievable  cancellation  ratio  as  a  function  of  the  power 
series  expansion  coefficients  for  the  two  nonlinearities,  {cx^}  and  {$.}.  Based 
on  Section  Al,  we  will  model  the  envelopes  at  the  n-th  harmonic,  v^  ntt)  and 
v2  n^)*  as  t*,e  sum  °f  their  principal  terms  plus  their  principal  corrupting 
terms. 


V  (t)  -  V  OL  +  -(-y2)-  a  V2(t) 
l,n  2n  n  4  n+2 


v  (t)  =  V  -&1  3  +  g,  V2(t) 

2,n  2n  n  4  n+2 


(A- 7) 


The  output  residue  power  from  a  cancellation  system  such  as  that 
shown  in  Figure  A-l  is 


I  e  I  2  =  / 1  V1(o))-WV2  (to)  |  2  dco 


•  where  V-^(u))  and  Vy(u))  are  the  Fourier  transforms  of  v-^(t)  and  vo(t)  and  W  is 
a  complex-valued  weight.  The  residue  |e|2  is  minimized  when  W  i  selected 
such  that 


giving 


/[V1(u))-WV2(ui)]V*(<jj)  du  =  0 

Jv^u^V^uOdai 


W  = 


/|v2(u))|2do) 

2 

Using  this,  |e|  becomes 

|e|2  »  /|v1(o))|/daj  - 


2  _  . . .  /IViWVjMl2^ 

/  I V 2  (to )  |  2dw 


Normalizing  this  by  the  power  in  the  unweighted  channel  yields  the  cancella¬ 
tion  ratio 


CR 


=  1 


/ 1  V1(o>)  |  2du) 

/ 1  V1<to)V*<oo)  |  2doj 


J|v1(u))!2dwJ!v2(w)|2du) 
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FIGURE  A-l .  MODEL  FOR  CANCELLATION  ANALYSIS 


Applying  Parseval'  Theorem  and  retaining  only  the  envelope  terms  yields  the 
optimum  cancellation  ratio 


CR  =  1  - 


[/v1(t)v2(t)dt]2 
Jv2(t)dt  Jv2(t)dt 


(A-8) 


where  CR  is  the  numerical  ratio  of  the  output  power  to  power  at  the  unweighted 
input.  Inserting  (A-7)  into  (A-8)  gives 


CR  =  1  - 


i/[avn(t)  + 


,vn+2(t)][e  vn(t>  +  n+2°  --n+2 


n  '  4  n+2  it,  +  ~^n+2 


vn+/(t)]dt}‘ 


+  ar<W'n4'2(0]2dt  /»„*”(«)  +I=|aen+2v"+2<')]2dt 
♦  iS±»taii8n+2+SnCln+2)v2n+2(t)  + 

^-2"+2<t,  +  <s«42  2n+4(t)]dt  -  -(A' 


(A- 9) 


/[ey"<t)+en6ii+2is±21v2»«(t)  +  isg42+2v2n+4(t)]dt 


As  a  worst  case  let  us  consider  the  AM  transmitter  to  be  sinusoi¬ 
dally  modulated,  double  sideband  suppressed  carrier.  This  is  chosen  as  a  worst 
case  because  the  envelope  covers  the  full  range  between  its  peak  value  and  zero, 
and  because  it  spends  a  large  percentage  of  its  time  near  its  peak  value 
where  the  harmonics  generated  are  strongest.  Thus,  we  write 


V(t)  =  A  cosu)  t 
m 


(A-10) 


The  evaluation  of  the  integrals  in  (9)  is  straightforward,  resulting 
in  the  following 

r„2 _  A2n  ,2n. 


Jv  n(t)dt  -  ^  O 


rv2n+2 ,  . ,  _  A2n+2  2n+2. 

jV  (t)dt  ^2  (  n+1 

r_2n+4.  ,,  A2n+<i  ,2n+4. 

jv  (t)dt  =  (  n+2) 


A  (2n+l) 
2 (n+1) 


.2n  0 

K-  <  n> 

22n  n 


A  (2n+l)(2n+3) 
4 (n+1) (n+2) 


.  2n  0 

V  <  ;>  (A’n) 

22n  n 


66 


(A- 12 


where  the  notation  (J)  is  the  binomial  coefficient 

(k)  -  - - 

Y  j ! (k-j ) ! 

Inserting  (11)  into  (9)  gives 

r  p  (n+2) rA2(2n+l), ,  „  ,«  >>  .  (n+2) 2 fA4 (2n+l) (2n+3) 

i0lnBn  +  4  2 (n+1)  J  (0lnBn+2+BnCtn+2)  +  4^+1)  (n+2  )' 

an+2Bn+2^ 


r  2  (n+2)  rA  (2n+l)  (n+2T2  rA^(2n+l)  (2n+3), ,  , 

lWn+2  2  1  2  (n+1)  J  +  16  an+2 l4 (n+1) (n+2) 

r„2.„  „  (n+2)  ,A2(2n+l),  .  (n+2)2„2  rA4(2n+l) (2n+3) 

^ n+BnBn+2  2  1  2  (n+1)  J  +  ~16“Bn+2  YTn+lT(n+2) 

Collecting  terms  in  (13)  over  a  common  denominator  results  in 

A4 (a  a  . „)2 (n+2) (2n+l) 


i  i  a  a  .  0A  (n+2)  (2n+l)  a  .,A  (n+2)  (2n+l)  (2n+3) 

64<"+1>  <“n  +  ---\(n+l) - +  — --rtlSn - }* 

,  8  8  ~A2 (n+2) (2n+l)  $2  A4 (n+2) (2n+l) (2n+3) 

lPn  +  4 (n+1)  64 (n+1)  ' 

Dividing  numerator  and  denominator  by  a%2  and  dropping  higher  order  terms 
results  in 

.4  Bn+2  an+2N2,  ,  „w„  ,  ,x 
A  -5 - - — )  (n+2)  (2n+l) 


-  x  2  r-.  .  n+2  A‘  (n+2 )  C2n+1 )  T .  r  n+2  A~ (n+2) (2n+l) , 

64 (n+1)  tl  +— - 4?S+1) - }  U  - 405+1)  1 


Note  that  if  the  nonlinearities  are  identical,  the  numerator,  and  hence  the 
CR,  go  to  zero,  implying  perfect  cancellation. 


ACHIEVABLE  CANCELLATION  RATIO  IN  TERMS  OF  MEASURABLE  HARMONIC 
POWER  LEVELS 


The  coefficients  {a^}  and  are  not  easily  measured  directly, 

and  hence  the  utility  of  Equation  (A-15)  is  less  than  ideal.  However,  these 
coefficients  may  be  directly  related  to  the  power  measurements  made  at  the 
individual  harmonics  with  the  same  AM  fundamental  to  provide  a  more  useful 

result. 


We  approximate  the  power  levels  at  the  n-th  harmonic  of  each  non¬ 
linearity,  by  using  only  the  first  term  of  (A-7)  to  obtain 


Pl»n  = 

/Vl2,n 

P2  ,p.  = 

/V2,n 

n.2  r„2n . 


n\2  r„2n, 


(A-16) 


Performing  the  integrations  with  V(t)  given  by  (A-10) ,  we  obtain 

2.2n 


>  ~ 

a  A 
n 

<2"> 

l,n 

2  An 

i  ~ 

e2A2n 

n 

<2:> 

2,n 

24n 

power 

levels 

at  the 

• 

2 

„  an+2 

A2(n+2) 

l,n+2 

2  A (n+2) 

„  0n+2 

A2(n+2) 

2  ,n+2 

2 A (n+2) 

<2"Th 


(A-17) 


(A-18) 


The  ratios  of  the  corresponding  expressions  in  (A-17)  and  (A-18)  is 


.  .2n+4. 

l.n+2  „  ,an+2,2  AJ  }  n+2}  _  ,an+2,2  A 4  (2n+l)(2n+3) 

P_  V  n  *  ft 


P,  '  a  '  08  ,2n.  v  a 

1  ,n  n  2  (  )  n 

n 

p  g  / 

2  ,n+2  _  /jn+2.2  £_  (2n+l)  (2n+3) 
P2,n  0n  ^  26  (n+D(n+2)  ’ 


6  (n+1) (n+2) 


(A-19) 
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whence 


I  o 

- - - 1  ‘/1Pl,n+2(n+1)(n+2)J/LPl,n(2n+1)(2n+3)J 

Bn+2A2  3  _ 

B  =  2  ^fP2,n+2^n+1^  ^n+2^]/{P,  (2n+l) (2n+3) ] 

We  now  insert 

(n+2)2(/p  /p  -  /p"  ~JV  )2 

_ 2,  n+2  2,n _ l,n+2  l,n _ _ 

(n+1)  (2n+3)(l  +  2*/[P;L  n+2/P1  n]  *  [  (n+2)3(2n+l) /(n+1)  (2n+3)]  }  • 


•  {l+2/[P2>n+2/P2jn]-[(n+2)3(2n+l)/(n+l)(2n+3)]} 
For  large  n,  we  can  further  simplify  (A-21)  to  obtain 


(A-20) 


(A-21) 


</P2.n+2^P2.ti  ~  /Pl,n+2^Pl.n* _ 

2(1+2”-/Pl,n+2/pi  ,„J<1+2"*2  ,„+2/P2,n)  ’ 


for  n-5. 


(A-22) 


3.1  EVALUATION  FOR  SOME  SPECIAL  CASES 

Equation  (22)  will  be  evaluated  for  two  cases. 

Case  1 :  Rapid  decrease  in  harmonic  energy. 

Suppose  2n^lin+2/Pl  n  <  1  and  2n^^7f^;  <1.  Then 

'  'Pl,n«/Pl,„  <  2'2" 

and 

as  an  upper  bound. 

Case  2 :  Slow  decrease  in  harmonic  energy. 

Suppose  2»/P1>lrt2/P1(Il  >  1  end  2n/p^7p^  >  1.  Then 

“  =  ^2<Al,n/Pl.«P2  '  ^2.n/P2,W2>2 

whence  we  can  upper-bound  CR  as 
CR  <  1/2. 
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V. 


These  results  lend  credence  to  the  intuitive  conclusion  that  where 
the  harmonic  power  levels  deer... use  slowly,  the  achievable  cancellation  ratio 
is  much  more  dependent  on  mulching  ui  the  noniir.earit ies .  Case  1  above  cor¬ 
responds  to  the  situation  where  the  harmonic  level  falls  off  rapidly.  Case  2 
corresponds  to  the  situation  where  the  harmonic  level  does  not  decay  as  rapidly, 
thereby  leading  to  a  smaller  cancellation  ratio  as  shown  by  the  upper  bound 
found  above. 


APPENDIX  B 

LISTING  OF  FORTRAN  PROGRAM 


This  appendix  contains  a  listing  of  the  FORTRAN  program 
simulates  the  key  portions  of  an  harmonic  cancellation  system. 


which 


f 


C  HARMONICS  CANCELLATION 

5017  FORMAT!  2X,’ Si.  1UL  PROG;  NLN01  FUND;  Ail/SC  ’) 

DIMENSION  At  10)  ,  1UW12100) 
l,IOi;(2100),ALHK(10),AU»<(  10) 

COMPLEX  CVO,CVu,CLRR,ClU,CCSU,CV 

c 

TI*2.5 
NO*  1  1 
FS= 10000. 

C=1 .0 
FO* 1217 , 

1  TYPE  501 1  ,Tl  ,NC,FS,0,FO 

5011  FORMAT (2X, '  RUM  TiMEC  1  )=  '  ,F10.4  , '  NUM  AUDR  >SITS(  2  )=  ’  ,  lA  , 

1  /  ,  2X  ,  '  S AMPLE  RATt.(  3)*'  ,F9. 0,  '  CAI.\(4  )='  ,  F  H).  3  , 

2  / , 2X  , '  FUND  FRE0(5)  =  *F9.0) 

TYPE  5013 

5013  FOR;  IAT ( 2X , *  ENTER  DESIRED  CHANGE :  RUM  AND  VALUE') 

ACCEPT  * ,  ICliKG  ,VaL 
1CENG= 1CHNC+1 

CO  TO  (70.10,15,20,25,30)  ICIi.YC 
10  T1=VAL 

GO  TO  1 
1 5  i\'G=VAL 

GO  TO  1 
20  FS*V'AL 

GO  TO  1 
25  G=VAL 

GO  TO  1 
30  FU»VAL 

GO  TO  1 

70  CALL  ASSIGN  (fa  , 9I1SAXNL.  OAT ,  9  >  ,  ’  CG'  ) 

5007  FORMAT  (2X,'SMPL  RATE  =  '  ,F9 . 0 , 3X  , '  NUM  AJ)OK  bITS= ’  , 14 , JX  , 

1  'GAIN-' ,F10.3) 

5009  FORMAT  ( 2X, '  FUND  FKE()*' ,F9.0,3X, ‘RUN  TIME* '  , F10.4 , 3X) 

TYPE  5017 

PRINT  5017 

PRINT  5007,  FS,NG,G 

PRINT  5009,  FO , T I 

DO  150  I A- 1,2 100 

IAW(IA)=0 

iqw(iA)»i) 

150  CONTINUE 
N»3 

Pl=3. 14159 
PI2»Pl/2. 

PI4=PI/4. 

P 14 j»3.*PI4 
ERS<)S*0.0 
HRSgS-0.0 
K*  1 

W0-2. *PI*KO 


»■  a»- 


JF-0 

DLL  1  - 1  /  FS 

t-dclt 

FLI M=2 . 0**NC- 1  .  U 
100  CONTI.N UK 

AVu=0.98*:;iN(00*i) 

ANGll=Pl4 

IKAVO.LT.O.)  AN  OH  =  -PI43 
CVu=CMPLX( AVO ,0.0) 
c 

AVIi*AV0**2 

AVtlu=l  +(AVIi/ 4.  )*(1  .  +  (AVli*.  1  )*(  1  .+AVI1/18.  )  ) 
AVIln«0. 77563*AV11ii*AVU*aVO 
AViID=0.98*aVO**5 
CVll*C.TPLX(AVlili,  AVlUi) 

C 

CV0=rLIM*CV0 
CVll*KLl!l*CVti 
AC VO=CAHS ( CVO ) 

KCCSI1=C0S(.\NCI1) 

qccsh=si:.(akgh) 

CCSli=l  .4l4*CMPLX(KCCSL,gCCSll) 

11T  ACVO.GC.  KLIM)  AC V0=  (  FLI.T) 

IA=IF1X(ACV0)+1 

KSii  =  FL0«\T(  IRW(  IA) )*UEAL( CCSU) 

QSH=FLOAT(IgW(IA))*AlMAC(CCSll) 

RERK*-KEAL(CVU)+USH 

QERK=-AI!tAC(CVU)+QSi; 

rcek=-g*rlrr*sign(i  .  .rccsr) 

QCLR»-G*QERR*SIGN(1.  .GCCSil) 

IF  I  RCER.GE .  FLIM)  RCE  R»  F  L I  f  1 
I F( KCER.  LL . -F'LI.'I)  RCLR--FLIM 
IFiQCEK. CE . FLI.'l)  qCEK=FLIM 
lF(gCi:K.LE.-FLIM)  gCEK=-FLlli 
igCEk* I F 1 X (QCEU ) 

Ii<CER=l  F  1X(  RCEK  ) 

I.*W(  IA)=IRW(  IA)+IKCEK 
igW(IA)»IQU(I.\)+lgCbR 
I F(  I (  1  A)  . GE .  FLIM)  IKW(  IA)-FLLI 
IF  (IKW(Ia).  LE  .-FLIT)  ll<W(  1A  )  —  FLU 

if(iqw(ia).cf:.fli  i)  igw(iA)-FLi;i 

IF (  Ii')U(  1A)  .LE.-FLIM)  IQW(  IA)  —  FLIM 

LKSg«KEKK**2+<)CKR**2 

ilARSQ-(Ri:AL(CVU))**2+(ALTAC(CVll))**2 

EllSQS-ERSQS+ERSq 

iiKSQS=UKSQS+ilARSg 

IF(K.Eg.lOO)  GO  TO  300 

K=K+ 1 

GO  To  250 
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300 


K-l 

aer=ersqs/hrsqs 

AtiH=URSQS/ 100.0 
EKSgS=0.0 
klkSQS-0.0 

1F(AER.LT.  IE-9)  AEK-1E-9 
1F(A11R.LT.  IE-9)  AHR-1K-9 
J  P “MOD ( J  P ,  1 0 )+ 1 
ALEK( JP)=10*ALuG10( \EU) 

ALHR(  JP)”10*ALOC10(AiiR) 

IF(JP.NE.  10)  Gi i  TO  250 
TYRE  222  ,ALKR,T  ALliK 
PRINT  222  ,ALEK,T,ALJ1R 
222  FORnAT(4X,10(F6.1,lX),2X,F7.5/4X,10(F6.1,lX)) 

250  CONTINUE 

255  T-DELT+T 

IF(T.CT.TI)  GO  TO  999 
GO  TO  100 
999  CONTINUE 

1  MAX— 14000 

1M IN- 14000 

DO  650  J-1,2100 

IMAX-tlAXO( IMAX ,  1RW( J)  ,  IQW(J)  ) 

650  1MIN-MIN0( IMIN , IRW( J ) , IQW( J ) ) 

MIMIN-IABS(IMIN) 

L'IAX»MAXO(MIMIN ,  IMAX) 

PRINT  5021, IMAX 
TYPE  5021,  IMAX 

5021  FORMAT ( 2 X,* MAX  MEM  VALUE  -  ’,110) 

TYPE  5000 

5000  FORMAT(2X, 'ENTER  1  FOR  COMPLEX  WEIGHT  PLOT...’ 
ACCEPT  5001 , IPLOT 

500 1  FORMAT  (II) 

1F( IPLOT. NE.l)  GO  TO  1100 
FSCLF-230.0/IMAX 
DO  675  KA-4,2100,4 
J-KA/4 

675  IRW(J)=IRW(KA)*FSCLF 

DO  690  L-4,2100,16 
JA-L/4 

IQW(JA)-IQW(L)*FSCLF 

LA-L+4 

IQW( JA+J )-IQW(LA)*FSCLF 
IQW(JA+2)-0 


t>90 


701) 


5003 


1100 


IQW(JA+3)-0 
CALL  PLOT55(9 ,0 ,0) 

CALL  PLOT55 ( 10 , ,  ) 

CALL  PLOT55( 2, 1+2+4+32+64+512,) 
CALL  PLOT55(7 ,0,0) 

CALL  PLOT55( 1,0,) 

1X-0 

IY-0 

DO  700  J-1,11 
CALL  PLOT55(5 , IX  ,  1 ) 

CALL  PLOT55 (4 . 1 , IY) 

IX-IX+50 

IY-1Y+23 

CONTINUE 

CALL  PLOT55(3 ,-512, 1KW ) 

CALL  PLOT55( 1 , 1 , ) 

CALL  PLOT55(3 ,-5 12 ,  1QW) 

READ(5 ,5003) 

FORMAT ( 15) 

CALL  PLOT55(2, 1+2+4+32+64+512,) 
CALL  PLOT55(2 ,512,) 

CALL  PLuT55(9 ,0 ,0 ) 

CALL  PLOT 55 (10,,) 

CALL  PLOT55(2,, 1+2+4+32+64+512) 
END 


75 


REFERENCES 


11]  "Nonlinear  Interference  Cancellation  System,"  RADC-TR-78-225 ,  November 
1978,  B032512L. 

12]  Higa,  Walter  H. ,  "Spurious  Signals  Generated  by  Electron  Tunneling  on 
Large  Reflector  Antennas,”  Proc  IEEE,  vol.  63,  no.  2,  February  1975. 

[3]  Chase,  W.M. ,  J.W.  Rockway  and  G.C.  Salisbury,  "A  Method  of  Detecting 
Significant  Sources  of  Intermodulation  Interference,"  IEEE  Trans  on 
Electromagnetic  Compatibility,  vol.  EMC-17,  no.  2,  May  1975. 

141  Chase,  W.M.,  "Ship  RFI  Survey  Procedure  for  HF  Frequencies,"  NELC 
Technical  Document  336,  21  June  1974. 

15 J  Wiener,  N.  Nonlinear  Problems  in  Random  Theory ,  MIT  Press,  1958. 

[ 6 ]  "Nonlinearity  System  Modelling  and  Analysis  with  Application  to 
Communica tion  Receivers,"  RADC  TR-73-178,  June  1978,  766278/6. 

]7]  Thomas,  E.J.,  "An  Adaptive  Echo  Canceller  In  a  Nonideal  Environment 
(Nonlinear  or  Time  Variant),"  BSTJ ,  vol.  50,  no.  8,  October  1971. 

18]  Falconer,  D.D. ,  "Adaptive  Equalization  of  Channel  Nonlinearities  in 
QAM  Data  Transmission  Systems,"  BSTJ,  vol.  57,  no.  7,  September 
1978. 

19]  Davenport,  W.B.,  Root,  W.L.,  An  Introduction  to  the  Theory  of  Random 
Signals  and  Noise,  McGraw-Hill  Book  Co.,  New  York,  1978,  p.  ~2847 

110]  Jolley,  B.W. ,  Summation  of  Series,  Dover  Press,  19bl. 

111]  M.  Schetzen,  "Determination  of  Optimum  Nonlinear  Systems  for 
Generalized  Error  Criteria  Based  on  the  Use  of  Gate  Functions,"  IEEE 
Trans  on  Information  Theory,  January  1965,  pp.  117-125. 


76 


4 


MISSION 
of 

Rome  Air  Development  Center 

RA VC  plank  and  execute*  research,  de.veZopme.nt,  tut  and 
selected  acquisition  programs  In  support  oi  Command,  Control 
Communications  and  Intelligence.  (C3I )  activities .  Technical 
and  engineering  support  within  areas  oi  technical  competence 
is  provided  to  ESP  Program  OHices  IPOs)  and  other  ESP 
elements.  The  principal  technical  mission  areas  are 
communications,  electromagnetic  guidance  and  control,  sur¬ 
veillance  oi  ground  and  aerospace  objects.  Intelligence  data 
collection  and  handling,  iniormatlon  system  technology. 
Ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 


